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Recently it has been shown that when an equation that allows so-called pulled fronts in the mean- 
field limit is modelled with a stochastic model with a finite number TV of particles per correlation 
volume, the convergence to the speed v* for TV — *• oo is extremely slow — going only as In -2 TV. 
Pulled fronts are fronts that propagate into an unstable state, and the asymptotic front speed is 
equal to the linear spreading speed v* of small linear perturbations about the unstable state. In 
this paper, we study the front propagation in a simple stochastic lattice model. A detailed analysis 
of the microscopic picture of the front dynamics shows that for the description of the far tip of the 
front, one has to abandon the idea of a uniformly translating front solution. The lattice and finite 
particle effects lead to a "stop-and-go" type dynamics at the far tip of the front, while the average 
front behind it "crosses over" to a uniformly translating solution. In this formulation, the effect 
of stochasticity on the asymptotic front speed is coded in the probability distribution of the times 
required for the advancement of the "foremost bin". We derive expressions of these probability 
distributions by matching the solution of the far tip with the uniformly translating solution behind. 
This matching includes various correlation effects in a mean-field type approximation. Our results 
for the probability distributions compare well to the results of stochastic numerical simulations. 
This approach also allows us to deal with much smaller values of TV than it is required to have 
the In -2 TV asymptotics to be valid. Furthermore, we show that if one insists on using a uniformly 
translating solution for the entire front ignoring its breakdown at the far tip, then one can obtain 
a simple expression for the corrections to the front speed for finite values of TV, in which various 
subdominant contributions have a clear physical interpretation. 

PACS Numbers: 05.10.Gg, 05.40.-a, 05.70.Ln, 82.20.Mj 



I. INTRODUCTION 
A. Fronts and fluctuation effects 

In pattern forming systems quite often situations oc- 
cur where patches of different bulk phases occur which 
are separated by fronts or interfaces. In such cases, the 
relevant dynamics is usually dominated by the dynamics 
of these fronts. When the interface separates two thermo- 
dynamically stable phases, as in crystal-melt interfacial 
growth problems, the width of the interfacial zone is usu- 
ally of atomic dimensions. For such systems, one often 
has to resort to a moving boundary description in which 
the boundary conditions at the interface arc determined 
phenomenologically or by microscopic considerations. A 
question that naturally arises for such interfaces is the 
influence of stochastic fluctuations on the motion and 
scaling properties of such interfaces. 

At the other extreme is a class of fronts that arise in 
systems that form patterns, and in which the occurrence 
of fronts or transition zones is fundamentally related to 
their nonequilibrium nature, as they do not connect two 
thermodynamic equilibrium phases which are separated 
by a first order phase transition. In such cases — for ex- 
ample, chemical fronts the temperature and density 
transition zones in thermal plumes , the domain walls 
separating domains of different orientation in in rotating 
Rayleigh-Benard convection |l , or streamer fronts in dis- 



charges jlj — the fronts are relatively wide and therefore 
described by the same continuum equations that describe 
nonequilibrium bulk patterns. The lore in nonequilib- 
rium pattern formation is that when the relevant length 
scales are large, (thermal) fluctuation effects are rela- 
tively small For this reason, the dynamics of many 
pattern forming systems can be understood in terms of 
the deterministic dynamics of the basic patterns and co- 
herent structures. For fronts, the first questions to study 
are therefore properties like existence and speed of prop- 
agation of the front solutions of the deterministic equa- 
tions, which in most cases are partial differential equa- 
tions. In the last decades, the fundamental propagation 
mechanism of such deterministic fronts has become rela- 
tively well understood. 

From the above perspective, it is maybe less of a sur- 
prise that the detailed questions concerning the stochas- 
tic properties of inherently nonequilibrium fronts have 
been addressed, to some extent, only relatively recently 
P-pl, and that it has taken a while for researchers to 
become fully aware of the fact that the so-called pulled 
fronts p|-pT[ which propagate into an unstable state, do 
not fit into the common mold: they have anomalous sen- 
sitivity to particle effects f§-|ll| , and have been argued 
to display uncommon scaling behavior [ p^ , |l9| -^l| . 

Pulled fronts are fronts which propagate into an unsta- 
ble state, and whose propagation dynamics is essentially 
that they are being "pulled along" by the growth and 
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spreading of the small perturbations about the unstable 
state into which the front propagates — their asymp- 
totic speed v as is equal to the linear spreading speed v* of 
perturbations about the unstable state: v as — v* jl4]-{l7|. 
This contrasts with the pushed fronts, for which v as > v* , 
and whose dynamics is determined by the nonlineari- 
ties in the dynamical equations [|il^-|l7||. The behavior 
of pushed fronts is essentially similar to fronts between 
two (meta)stable states. 

The concept of a pulled front most naturally fits a for- 
mulation of the dynamical equations in terms of con- 
tinuum variables, for by "small perturbations" we mean 
that the deviations of the field values from the values 
in the unstable state are small enough that nonlinear 
terms in the deviations can be neglected. From vari- 
ous directions, it has become clear in the last few years 
that such pulled fronts do show very unusual behavior 
and response to perturbations. First of all, Brunet and 
Derrida have shown that when the continuum field equa- 
tions are used for a finite particle model so as to have 
a growth cutoff at the field value 1/N, where N is the 
typical number of particles in the bulk phase behind the 
front, the deviation from the continuum value v* of the 
front speed is often large, and it vanishes only as 1/ In TV 
(with a known prefactor which they calculated) ||. On 
the other hand, we recently found that with an infinites- 
imal growth cutoff and a similarly infinitesimal growth 
enhancement behind it, one can have a much higher front 
speed than v* Jig] . Furthermore, the scaling properties 
of pulled fronts in stochastic field equations with a par- 
ticular type of multiplicative noise have been found to 
be anomalous: in one dimension, they are predicted to 
exhibit subdiffusive wandering JlSj ] , but in higher dimen- 
sions their scaling behavior is given by the KPZ equation 
p2| in one dimension higher than one would naively ex- 
pect [^0|,|lj (the question to what extent these results 
are applicable to lattice models, where the finite particle 
effects always make the fronts weakly pushed, is still a 
matter of debate (HJ2J]). Moreover, even without fluc- 
tuations, pulled fronts respond differently to coupling to 
other fields, e.g., they never reduce to standard moving 
boundary problems, even if they are thin p5[ . 

All these effects have one origin in common, namely 
the fact that the dynamics of pulled fronts, by its very 
nature, is not determined by the nonlinear front region 
itself, but by the region at the leading edge of the front, 
where deviations from the unstable state are small. To 
a large degree, this semi-infinite region alone determines 
the universal relaxation of the speed of a deterministic 
pulled front to its asymptotic value ||,[l(|[l7]] , as well 
as the anomalous scaling behavior of stochastic fronts 
fL9| ^l| , p3 24 1 in continuum equations with multiplicative 
noise. As realized by Brunet and Derrida ||, the crucial 
importance of the region, where the deviations from the 
unstable state are small, also implies that if one builds a 
lattice model version of a front propagating into an un- 
stable state, the front speed is surprisingly sensitive to 
the dynamics of the tip (the far end) of the front where 



only one or a few particles per lattice site are present. It 
is this effect which is the main subject of this paper. 



B. Open questions 

If we study fronts for a field describing the number 
density <f> of particles, and normalize the field in such a 
way that its average value behind the front, where there 
are N particles per unit of length, is 1, then at the very 
far end of the leading edge, where the discrete particle 
nature of the actual model becomes most noticeable, the 
value of the normalized number density field is of order 
1/N. Brunet and Derrida |j| therefore modeled the effect 
of the particle cutoff in their lattice model by studying 
a deterministic continuum front equation, in which the 
growth term was set to zero for values of <f> less than 
1/N. They showed that this led to a correction to the 
asymptotic front speed of the order of 1/ In 2 N with a 
prefactor, which is given in terms of the linear growth 
properties of the equation without a cutoff. Because of 
the logarithmic term, in the dominant order, it does not 
matter whether the actual cutoff should really be exactly 
1/7V (corresponding to exactly one particle), or whether 
the growth is just suppressed at values of <f> of order 1/N, 
since 1/ In 2 (cN) w 1/ln 2 A in dominant order. Simula- 
tions of two different lattice models by Brunet and Der- 
rida H and by van Zon et al. [^6| gave strong support for 
the essential correctness of this procedure for sufficiently 
large N, but showed that there can be significant devi- 
ations from the asymptotic result for large but not ex- 
tremely large N. Moreover, for a different lattice model, 
Kessler et al. ]l0| did observe a correction to the average 
front speed of order 1 / In 2 N but with a prefactor which 
they claimed was a factor of order two different from the 
prediction of Brunet and Derrida. 

There are hence several questions that lead us to re- 
consider the finite particle effects on the average front 
speed of pulled stochastic fronts: 

(i) Why is it that a simple cutoff of order 1/N in 
a deterministic equation for a continuum (mean- 
field type) equation apparently leads to the proper 
asymptotic correction to the average speed of a 
stochastic front? 

(ii) Can we get a more microscopic picture of the 
stochastic behavior at the far end of the front, 
where there are only a few particles per lattice site? 

(in) Can we go beyond the large A asymptotic result 
of Brunet and Derrida, e.g., can we calculate the 
correction term for large but not extremely large 
values of N or even for arbitrary A? — After all, 
one might a priori expect correlation effects to be 
very important for fronts whose propagation speed 
is strongly affected by the region where there are 
only a few particles per site. 
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(iv) What is the role of correlation effects? 

(v) To what extent do the specific details of the partic- 
ular stochastic model play a role? 

(vi) Can one resolve the discrepancy noted by Kessler 
et al. 0? 



C. Summary of the main results 



In this paper, we address these questions and answer 
the majority of them for a specific model for which 
Breuer et al. || already studied the asymptotic speeds 
of stochastic fronts numerically a number of years ago. 
The model consists of particles making diffusive hops on 
a one-dimensional lattice, and being subject to growth 
and death on each lattice site. It is very close to the 
one also studied by Kessler et al. p0| , the only differ- 
ence being that their model includes a correlation term, 
which is small and irrelevant for large N. The absence of 
such correlations makes the model studied by Breuer et 
al. easier to analyze. Moreover, an examination of the 
numerical results therein shows that the deviation of the 
asymptotic front speed from its pulled front value indeed 
behaves as 1/ In 2 N, (although it was not realized in ||), 
with a prefactor that, over the range of TV-values stud- 
ied, is different from the one predicted later by Brunct 
and Derrida |J, but not as much different as Kessler et 
al. claimed it to be for their own model JlO[] . For each 
stochastic realization of a front, which moves into a re- 
gion where no particles are present, one can always iden- 
tify a foremost occupied lattice site or "the foremost bin" . 
In the region near this one, fluctuations are large and the 
discreteness of the lattice and of the particle number oc- 
cupation is extremely important: the standard descrip- 
tion, which assumes that the average particle density is 
uniformly translating, breaks down in this region. More- 
over, since the particle occupation numbers are small in 
the tip, essentially all known methods fail, based as they 
are on large- TV expansions. 

With a finite particle cutoff, fronts are never really 
pulled, but instead are weakly pushed p7j . Even for the 
simplest case of a pushed front in a second order non- 
linear partial differential equation, in general, the speed 
cannot be calculated explicitly. It should therefore come 
as no surprise that with the various additional complica- 
tions described above, we do not have a full first princi- 
ples theory that gives the front speed for finite values of 
N for the model we study. However, in this paper we do 
obtain a number of new results for the behavior in the 
far tip of the front as well as for the effect of the region 
behind the tip on the finite-iV corrections. These results 
can be tested independently and our numerical simula- 
tions largely support the picture that emerges from this 
approach. In terms of short answers to the questions 
raised above in Sec. IB, we find that 



(i) For extremely large N, the asymptotic results of 
Brunet and Derrida based on a simple cutoff of 
order 1/N in a deterministic equation for a con- 
tinuum (mean-field type) equation become essen- 
tially correct because all the essential changes are 
all limited to a few bins behind the foremost one, 
where the particle numbers arc finite and small; to- 
gether with the fact that 1/ ln 2 (ciV) w 1/ ln 2 (iV) to 
dominant order, this ensures the correctness of the 
asymptotic expression for N — > oo. 

(ii) Yes, one can get a more microscopic picture of what 
happens near the foremost bin of the front; we de- 
velop mean-field type expressions for the probabil- 
ity distribution that describes the "stop-and-go" 
type behavior there (Sec. IV), and show that the 
results compare well with numerical simulation re- 
sults Sec. V. 

(Hi) A first-principles theory for the stochastic front 
speed for arbitrary N seems virtually impossible, 
except possibly in some special limits, as in prin- 
ciple, it will involve matching the approximately 
uniformly translating average profile behind the tip 
of the front to the non-uniformly translating near 
the foremost bin, where standard methods do not 
seem to apply. 

(iv) Correlation effects are very important near the tip; 
we identify two of them and model one: rapid 
succ essive forward hops of the foremost particle, 
Sec. [VC1 and j umping back of the foremost par- 
ticle , Sec7|VC|. 

(v) The details of the particular stochastic model 
play a role for the corrections in the asymptotic 
front speed through the global average front profile 
(quantified by A of Sees. Ill and IV) and through 
the effective profile near the tip, but their effects 
are truly minute. We demonstrate this by means 
of a mean-field theory that tries to extend the uni- 
formly translating front solution all the way to the 
far tip of the front (described in Sec. 1IIC). In this 
theory, there is a quantity a associated with the ef- 
fective profile at the tip, and we show that these 
two quantities, A and a, provide only subdominant 
corrections to the asymptotic large N result. 

(vi) The model considered by Kessler et al. [|l0| is 
slightly different from the one considered by Breuer 
et al., in the sense that number of particles of each 
species is finite. However, a priori, one expects that 
this difference in the two models would not affect 
the speed corrections for large iV. Our own simula- 
tions confirm this, and give no sign of a discrepancy 
between the asymptotic large- N speed corrections 
obtained from the two models (Sec. VI). 



We finally note that in this paper, we will focus on the 
case where the growth and hopping terms for a few par- 
ticles are the same as those for a small but finite density 
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of particles. In such cases, the front speed converges for 
N — > oo to the pulled front speed v* of the corresponding 
mean-field equation. As we will discuss elsewhere [ jl8| , 
with only slight modifications of the stochastic rules for 
few particles, one can also arrive at situations where the 
limits do not commute, i.e., where the stochastic front 
speed converges to a speed larger then v* as N — ► oo, 
even though the stochastic model would converge to a 
the mean- field equation with pulled fronts in this limit. 



D. Complications Associated with Discreteness of 
the Lattice and Particle Numbers 

The challenge of understanding the propagation of any 
one of these fronts lies in the fact that as a consequence 
of the discrete nature of the particle events and of the 
particle number realizations, the natural description of 
the far tip is not in terms of a uniformly translating so- 
lution for the average number of particles in the bins (we 
call each lattice site a "bin"), but is in terms of discrete 
notions like the foremost bin, individual jumps, etc. An 
additional complication is that in the presence of fluctua- 
tions, the front position exhibits diffusion-like wandering 
behaviour, which have to be taken out in order to study 
the intrinsic stochastic front dynamics, just like capillary 
waves beset analy zing the intrinsic structure of a fluid in- 
terface (Sec. Ill B| ) . The implication of all this is that (i) 
in the presence of an underlying lattice, instead of being 
uniformly translating, the position of the foremost occu- 
pied bin advances in a discrete manner, and (ii) due to 
the discrete nature of the constituent particles, the posi- 
tion of the foremost bin advances probabilistically, as its 
movement is controlled by diffusion. 

Based on these ingredients and observations, the cen- 
tral theme in this paper revolves around a picture of the 
tip of the front that is totally different from the conven- 
tional picture of a pulled front. We present the picture 
here in terms of its simplified essence, as it is helpful for 
the reader to bear it in mind throughout this paper: we 
call the foremost occupied lattice site at the far end of 
the tip of the front "the foremost bin" . Therefore, the 
very definition of the foremost bin on a lattice site means 
that it is occupied by at least one particle and that all 
the lattice sites on the right of it are empty. Naturally, 
an empty lattice site (all the lattice sites on the right of 
which are also empty) attains the status of the foremost 
bin as soon as one particle hops into it from the left. In 
reference to the lattice, the position of the foremost bin 
remains fixed at this site for some time, i.e., after its cre- 
ation, a foremost bin remains the foremost bin for some 
time. During this time, however, the number of particles 
in and behind the foremost bin continues to grow. As 
the number of particles grows in the foremost bin, the 
chance of one of them making a diffusive hop on to the 
right also increases. At some instant, a particle from the 
foremost bin hops over to the right: as a result of this 
hop, the position of the foremost bin advances by one unit 



on the lattice, or, viewed from another angle, a new fore- 
most bin is created which is one lattice distance away on 
the right of the previous one. Microscopically, the selec- 
tion process for the length of the time span between two 
consecutive foremost bin creations is stochastic, and the 
inverse of the long time average of this time span defines 
the front speed. Simultaneously, the amount of growth 
of particle numbers in and behind the foremost bin itself 
depends on the time span between two consecutive fore- 
most bin creations (the longer the time span, the longer 
the amount of growth). As a consequence, on average, 
the selection mechanism for the length of the time span 
between two consecutive foremost bin creations, which 
determines the asymptotic front speed, is nonlinear. 

This inherent nonlinearity makes the prediction of the 
asymptotic front speed difficult. One might recall the dif- 
ficulties associated with the prediction of pushed fronts 
due to nonlinear terms in this context, although the na- 
ture of the nonlinearitics in these two cases is completely 
different. In the case of pushed fronts, the asymptotic 
front speed is determined by the mean-field dynamics of 
the fronts, and the nonlinearties originate from the non- 
linear growth terms in the partial differential equations 
that descr ibe the mean-field dynamics (As we discuss in 
Sec. [II B , if one does not take out the wandering of the 



front positions, then the nonlinear growth terms actu- 
ally do affect the stochastic front dynamics in a subtle 
way too). On the other hand, for fronts consisting of 
discrete particles on a discrete lattice, the corresponding 
mean-field growth terms are linear, but since the asymp- 
totic front speed is determined from the probability dis- 
tribution of the time span between two consecutive fore- 
most bin creations, on average, it is the relation between 
this probability distribution and the effect of the linear 
growth terms that the nonlinearities stem from. 

Our approach is to develop a separate probabilistic the- 
ory for the hops to create the new foremost bins, and 
then to show that by matching the description of the 
behaviour in this region to the more standard one (of 
growth and roughly speaking, uniform translation) be- 
hind it, one obtains a consistent and more complete de- 
scription of the stochastic and discreteness effects on the 
front propagation. In the simplest approximation, the 
theory provides a very good fit to the data, but our ap- 
proach can be systematically improved by incorporating 
the effect of fluctuations as well. Besides providing in- 
sight into how a stochastic front propagates at the far 
tip of the leading edge, our analysis naturally leads to 
a more complete description that allows one to interpret 
(though not predict) the finite N corrections to the front 
speed for much smaller values of N than that are neces- 
sary to see the asymptotic result of Brunet and Derrida 
H . As one might expect, for values of N where deviations 
from this asymptotic result are important, model-specific 
effects do play a role. 

For the major part of our analysis, we focus on the 
most relevant and illuminating case in which the diffu- 
sion and growth rates of the model are both of the same 
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order. This regime is the most illustrative as it displays 
all the aspects of finite particle and lattice effects most 
clearly. We also investigate the case when the diffusion 
rate is much smaller than the growth rate to illustrate 
the correlation effects. For all of these cases, the match- 
ing between the behaviour of the tip of the front and the 
standard description of a uniformly translating solution 
behind it is a complicated process, for the lack of a proper 
small parameter that allows one to do perturbation the- 
ory. 

The paper is organized in the following manner: in 
Sec. II, we describe our model (which is the same as in 
Ref. ||) and define the dynamics of the front. The crux 
of the paper is presented in Sec. IV, where we present a 
detailed analysis of the microscopic picture of the front 
dynamics and show that for the description of the far 
tip of the front, one has to abandon the idea of a uni- 
formly translating front solution. The lattice and finite 
particle effects lead to a "stop-and-go" type dynamics at 
the far tip of the front, while the average front behind it 
"crosses over" to a uniformly translating solution. In this 
formulation, the effect of stochasticity on the asymptotic 
front speed is coded in the probability distribution of the 
times required for the advancement of the foremost bin. 
We derive expressions of these probability distributions 
by matching the solution of the far tip with the uniformly 
translating solution behind. This matching includes var- 
ious correlation effects in a mean-field type approxima- 
tion. In Sec. V, we compare our theoretical predictions 
of Sec. IV with the stochastic simulation results. In addi- 
tion to that, in Sec. Ill, we argue that the corresponding 
front solution is a case of a weakly pushed front and ana- 
lyze an effective mean-field solution that extends all the 
way to the foremost bin (thereby ignoring its breakdown 
near the foremost bin). This allows us to rederive the 
asymptotic velocity expression of Brunet and Derrida |^] 
and obtain the further subdominant finite- N corrections 
to it. In Sec. VI, we carry out the full stochastic sim- 
ulation for the model considered by Kessler et ai, and 
finally, we conclude the paper with a discussion and out- 
look in Sec. VII. 



II. THE REACTION-DIFFUSION PROCESS X+Y 
^ 2X ON A LATTICE 

We consider the following reaction-diffusion process 
X+Y f± 2X on a lattice in the following formulation: 
at each lattice position, there exists a bin. We label the 
bins by their serial indices k, k = 1,2,3,. . . , M, placed 
from left to right. Each bin has an infinite supply of 
Y particles. An X particle in the fc-th bin can undergo 
three basic processes: (i) diffusion to the (k — I)-th or 
the (k + I)-th bin with a rate of diffusion 7. If an X 
particle in bin 1 jumps towards the left, or an X particle 
in the M-th bin jumps to the right, then they are im- 
mediately replaced, (ii) forward reaction to produce an 
extra X particle having annihilated a Y particle (X + Y 
— > 2X), with a rate 7 g , and (Hi) if there are at least two 
X particles present in the fc-th bin, then any two of the 
X particles can react with each other and annihilate one 
X particle to produce a Y particle (2X — > X + Y), with 
a reaction rate 7^. A state of the system at time t is 
given by the numbers of X particles in the bins, denoted 
as {Nx,N 2 ,...,N M ;t}. 

In the context of front propagation, the above model 
was first studied by Breuer et al §. Up to Sec. V of 
this paper, we will confine ourselves to this model only. 
In Sec. VI, we will consider a slightly modified version 
of this model, numerically studied by Kessler and co- 
authors (!(]] , in which the number of Y particles in any 
bin is finite, and the Y-particles can diffuse from any bin 
to its nearest neighbour bins with the same diffusion rate 

7- 

A. The Master Equation 

The discrete, microscopic description of the above 
reaction-diffusion process inherently introduces fluctua- 
tions in the number of X particles present in any partic- 
ular bin. This necessitates a suitable multivariate prob- 
abilistic description of the system. Let us denote the 
probability of a certain configuration {Ni, N 2 , . . . , Nm; t} 
at time t by P(Ni,N 2 , . . . ,NM',t)- The dynamics of 
P(N±, N 2 , ■ ■ ■ , Nm', t) is given by the following master 
equation: 



I 

N 2 , ...,N M ; t) = J2hK N k+i + l)P(Ni, N 2 , ...,N k - 1, N k + 1 + 1, N M ; t) 

k 

+ (V fe _ ! + 1) P(Ni, N 2 , . . . , Nk- 1 + 1, Nk — I, . . . , N M ; t) 
-2N k P(N 1 , N 2 , JV fc _i, N k , N k + 1 ,...,N M ; t) ] 
+ l g [(N k - l)P(JVi, N 2 , JV fc _i, N k - 1, N k + 1 , ...,N M ; t) 

-N k P{N u N 2 , N k - lt N k , N k + 1 , ...,N M ; t)] 
+ y [N k (N k + l)P(JVi,JV a ,... > JV fc _i,JV fc + l,N k+u ...,N M ;t) 

-N k (N k - l)P(Nx, N 2 , iVfc_i, N k , N k + U ...,N M ; t)}} . (2.1) 
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The above equation is actually not quite accurate at the 
1st and M-th boundary bins, but we refrain from writ- 
ing out the correction terms explicitly, as they are not 
needed in the analysis below. 



B. The Macroscopic Density Field and the 
Fisher-Kolmogorov Equation 

If the forward reaction rate, j g> is much larger than 
the annihilation rate 7^, an initial conglomeration of X 
particles will start to grow in size as well as in numbers. 
To study this growth phenomena, we define (Nk(t)), the 
average number of X particles in the fc-th bin at time t, 
as 



(N k {t))= N k P(N u N 2 , 

{-Nfc/} h /=i...jv 



N M ; t) . (2.2) 



Using Eq. (2.1), it is easy to obtain the time dynamics 
of (N k (t)}, given by 



0_ 

at 



(N k (t)) = 7 [(N k + 1 (t)) + (JV fc _!(*)) -2(N k (t))} 
+ 7 9 (N k (t)) - ? [(N 2 k (t)) ~ (N k (t)}} > (2-3) 



with 



(N 2 k (t))= J2 NiP(N l7 N 2 , 



N M ;t). (2.4) 



For the sake of simplicity, we define 7 =j/j g , f = j g t 
and N — 2y n /~fd, and reduce the number of parameters 
in Eq. (|]§, to have f| 







oj{N k (t')) = 7 [{N k+ + (jv*-i(0) - 2 W)) : 



+ {N k (t')) - [{N 2 k (t')) - (N k (t>))] . (2.5) 



Following the procedure in Ref. 

i [(Ni(t)) - (N k (t))] term in Eq. Q by ^(N k (t)) 2 



if one replaces the 
1 



and further defines a mean "concentration field" on the 
fc-th bin by int roducing the variable 4> k = (N k ) /N, then 
from Eq. ( |2.5| ), one arrives at the following difference- 
differential version of the Fisher-Kolmogorov equation for 
the reaction-diffusion process X+Y f± 2X on a lattice, 
given by [§ 



a_ 

dt 



Mt) = 7 [h+i{t) + <t>k-i{t) - 2cf> k {t)} 



+ Mt) - 4>l(t). (2.6) 



The original Fisher-Kolmogorov or F-KPP equation 
p8| , p9[ is a partial differential equation in continuous 
space and time. Notice that in these variables, the prop- 
erties of the propagating front depend only on two pa- 
rameters, N and 7. 



III. MEAN-FIELD APPROXIMATIONS TO THE 
PROPAGATING FRONT SOLUTION 



As mentioned earlier, in this section we do not consider 
the proper stop-and-go type dynamics of t he t ip; instead, 
as a continuation of mean-field equation (2.6) above, we 
describe the entire front by the uniformly translating pro- 
file. We then make a number of general observations 
concerning the uniformly translating front solutions in 
mean-field type equations for the average profile, from 
the perspective of the questions raised in the introduc- 
tion. A central result of the discussion will be an expres- 
sion for the finite- iV value of the velocity, which allows 
us to interpret deviations from the asymptotic results of 
H in terms of physical properties of stochastic fronts. 



A. Front propagation in the dynamical equation for 

4>k(t) 



From the point of v iew of average number of X particles 
in the bins, Eq. (|1|) has two stationary states. One of 
them, for which (N k ) — N V k, is stable. The other, for 
which (Nit) = Vfc, is unstable. This means that any 
perturbation around the unstable state grows in time un- 
til it saturates at the stable state value. In particular, if in 
a certain configuration of the system, the stable and the 
unstable regions coexist, i.e., (N k ) = N V k < k and 
(Nk) =0 Vfc > ki, with k\ > ko, then the stable region 
invades the unstable region and propagates into it. In 
other words, in due course of time, the boundary between 
these two regions, having a finite width, moves further 
and further inside the unstable region. For a wide range 
of initial conditions, the speed, with which this bound- 
ary moves into the unstable region, approaches a fixed 
asymptotic value, v as - Simultaneously, the shape of this 
boundary between the two regions, determined by the 
average number of X particles, (N k ), plotted against the 
corresponding bin indices k, also approaches an asymp- 
totic shape. This asymptotic shape, therefore, becomes a 
function of (k — v as t) for long times, and this well-known 
phenomenon is known as the f ront propagation. In the 
present context, Eqs. (2.5-2.6) provides us with an ex- 
ample of front propagation into unstable states. We will 
follow the usual convention that the front propagates to 
the right in the direction of increasing bin numbers. 

In the mean-field approximation (2.6), the average 



particle density field 4> k (t) obeys a difference-differential 
equation. The asymptotic speed selection mechanism for 
propagating fronts into unstable states has been a well- 
understood phenomenon for a number of years, and it 
has been realized by various authors j^-[ll],[l7| that the 
calculation of the asymptotic front speed on a lattice 
for the type of Eqs. ( [2~5| - |2.6| ) proceeds along similar 
lines as it does for partial differential equations. It is 
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well-known that for Eqs. fl2.5[]2.6| ), the selection mecha- 
nism for v as depends entirely on the region, where the 
nonlinear saturation terms ( [(/V|(t)) — (Nk(t))] /N or 
4>k(t)) are much smaller in magnitude than the corre- 
sponding linear growth terms [(Nk(i)) or 4>k(t)], i.e., the 
leading edge of the front, where the value of (f>k(t) is 
very small, such that <^(i) <C <pk(t). In this region, 
the nonlinear terms can be neglected, and after having 
used 4>k(t) = (f>(k - v as t) = <p(£), where £ = k - v as t is 



the comoving coordinate, Eq. ( |2.q ) reduces to a linear 
difference-differential equation, given by 

- = 7 [HZ + i) + <HZ - !) - 2#£)] 



HO- 



(3.1) 



If one neglects the fact that the microscopic X particles 
are discrete and assumes that <f>{£) goes to zero continu- 
ously for £ — > oo, then a natural candidate for the solu- 
tion o f (/>( £) in the linear difference-differential equation, 
Eq. (3.1) above, is <f>(£) = A exp[— z£], where z is a 
real and positive quantity. With this solution of <f>{£) in 
the so-called leading edge of the front, one arrives at the 
dispersion relation 



Vas (z) = 



27[cosh(z) - 1] + 1 



(3.2) 



Like the other ex amp les of fronts propagating into un- 
stable states, Eq. (3.2) allows an uncountably infinite 
number of asymptotic velocities depending on the se- 
lected value of the continuous parameter z. However, for 
a steep enough initial condition that decays faster than 
exp(— A£) in £ for any A > Zq determined below (hence, 
a unit step function obeys this condition), the observed 
asymptotic speed equals the so-called linear spreading 
speed v* , given by v* = v*(z ), where z is the value of 
z, for which the dispersion relation v as (z) vs. z has a 
minimum. 

The fact that v* defined in this way is nothing but the 
linear spreading speed, i.e., the spreading speed of small 
perturbat ions whose dynamics is given by the linearized 
equation (3.2), follows from a saddle point analysis of the 
asymptotic be hav ior of the Green's function for the lin- 
ear equation (3.2), see e.g. |17|| . The name pulled fronts 
stems from the fact that this linear spreading almost lit- 
erally "pulls" the nonlinear front with it, the nonlinear 
terms just giving rise to saturation behind the front. 



B. The deceptive subtlety of the mean-field 
approximation 

As we discussed above, in the pulled front regime, the 
front speed of a given problem is determined completely 
by the linear term in the dynamical equation. In going 
from the exa ct equation ( p.5[ ) to the mean field approx- 
imation (|2.6|), we appear, at first sight, to have ignored 



only a term linear in <fik of o rder 1/N [the second term 
between square brackets in ( |2.5| )1. Hence, naively one 
might expect the front speed to converge as 1/N to the 
asymptotic value v*(zq). We already know from the work 
of Brunet and Derrida however, that the convergence 
is much slower, namely as l/ln 2 iV. How can the two 
results be reconciled? 

The resolution of the paradox lies in the fact that in 
the mean field approximation we completely ignore the 
diffusive wandering of fronts. If we follow the evolu- 
tion of an ens emb le of fronts, their positions [defined, 
e.g., by Eq. (1.1) below] will fluctuate: the root mean 
square wandering of the fronts grows as \ft as for any 
one-dimensional random walker [p|Jl2|] . This means that 
in reality the ensemble average (N/Jt)) does not acquire 
a fixed shape in the frame moving with the average speed. 
Instead, the average profile (iVfc(t)) continues to broaden 
in time, although the front shapes for the individual real- 
izations reach an asymptotic shape (see Fig. 5 of Ref. jf| 
for an illustration). This has a severe consequence: we 
cannot simply assume that the (Nj*(t)) term is small in 
the leading edge of the profile where (JVfc(i)) is small, and 
replace it by (Nk(t)} 2 — the few members of the ensem- 
ble, which are relatively further ahead, do give significant 
contributions through this term in regions where (Nk(t)} 
is small. Thus, while Eq. (2.5) is exact and contains the 
fluctuation effects due to the root mean squar e wa ndering 
of the front, the mean- field approximation ( |2.6| ) throws 
out such effects completely. 

If, on the other hand, we look at the shape of a partic- 
ular front realization in the appropriate position, so that 
the front wandering is taken out, the mean- field equation 
does yield a reasonably good description of this (condi- 
tionally averaged) front profile in the range where the 
particle occupation numbers are large and (hence) where 
fluctuation effects are small. Additional information is 
needed, however, to calculate the front speed. 

In passing, we note that the situation is somewhat 
similar to the theory of fluid interfaces: capillary wave 
fluctuations wash out the average interface profile com- 
pletely, but on scales of the order of the capillary length, 
the mean field theory for the so-called intrinsic interface 
profile works quite well. 



C. The front speed correction for large TV 

The above observations already allow us to arrive at 
and extend the results of Brunet and Derrida Q from 
a slightly different angle than in their original work as 
follows. First of all, from the discussion above, we notice 
that even though a mean field approximation (2.6) does 
not work for the ensemble-averaged front profile, but for 
a given stochastic front realization, the mean field theory 
does apply to a good approximation in the bins, where 
the number of X particles are relatively large. These are 
essentially the bins that arc sufficiently behind the fore- 
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most bin, the rightmost bin in given stochastic realiza- 
tion, on the right of which all bins are completely empty. 
Nevertheless, as mentioned in the beginning paragraph 
of this section, we assu me t hat the uniformly translating 
front solution of Eq. (2.£) holds for the description of 
the front profile all the way up to the foremost bin for a 
given rea lization. Secondly, the actual front solution of 
Eq. (2.6) is a case of a weakly pushed front as opposed 
to being a truly pulled front |2?],|3(|. This can be under- 
stood in the following manner: notice that in any bin the 
forward reaction X + Y — > 2X does not proceed unless 
there is at least one X particle in that bin to start with. 
As for any given realization of the stochastic front, the 
front propagation on a lattice is tantamount to the dis- 
crete forward movement of the foremost bin by units of 
1 (which can happen only through the diffusion of an X 
particle from the foremost bin towards the righ t), in the 
uniformly translating front solution of Eq. ( |2.6| ) , the dy- 
namics of the tip of the front is diffusion dominated. This 
makes any given realization of the front weakly pushed as 



opposed to being truly pulled, and moreover, the asymp- 
totic speed Wat is expected to be < v* for a finite N. This 
indicates that if we want to build all these in the same 
frame as in the velocity selection mechanism for a pulled 
front, one has to allow complex values of the parameter 
z (see Eq. (^2) and the discussion thereabove). Further- 
more, the existence of a foremost bin requires that the 
front profile must have a zero a bin ahead of the foremost 
bin. Having combined all these together, and without any 
loss of generality, we now require that the front profile in 
the linear region of Eq. (2.6) is given by [p|- pd| , p^T7| ] for 
</>(£) for v N < v* 



</>(£) = Asm[zi£, + f3] exp(-z r £) 



(3.3) 



such that </>(£) has a node at the coordinate of the bin 
just ahead of the foremost bin. In Appendix |A|, we show 
how Eq. (^) can be used to determine the complex de- 
cay rate z in terms of N and other parameters, and from 
that we obtain the deviation of the front speed vn from 
V*. The front speed vn is given by 



d 2 v n 



I'N 



dz 2 



O(zf) 



dz 2 



InN 



, A 

- z + In — 
a 



In < sin 



7TZ 



LaJV+1 



(3.4) 



where, according to Eq. ( p.2[) 



dz 2 



7 cosh zq 



zo 



(3.5) 



In the limit of large TV, the above result (3.4) reduces to 



d 2 V as 



dz 2 



In 2 N 



(3.6) 



which is nothing but the asymptotic expression for the 
velocity correction derived by Brunet and Derrida 
Their approach is based on the partial differential equa- 
tion analog of the mean field dynamical equation (2.6), 
in they introduced an artificial cutoff for the growth term 
for values of </>(£) < e, where e ~ 1 /N, to mimic the dom- 
inant role played by diffusion at the tip of the front as 
opposed to the growth term. 



D. Implications and discussion 

The above expressions for the speed corrections are al- 
ready quite instructive. First of all, as w e pointed out, 
for the speed differ enc e v* — Vn, Eq. ( |3.5| ) reduces to the 
expression of Eq. ( |3.6| ) of Brunet and Derrida ||] at the 
dominant order in the limit of very large N. To this or- 
der, the speed change is given explicitly in terms of N. 
The more general expression, Eq. (3.5), however, con- 
tains the factors A and a; these affect the sMfedominant 



behaviour, i.e., the corrections to the asymptotic large 
N expression. For realistic values of N, the corrections 
to the asymptotic behavi or ca n be quite significant 
As we shall show in Sec. 



VB, A depends on the global 



behaviour of the average front solution, including the be- 
haviour in the region where nonlinearities are important. 
This makes its value vary from model to model and it is 
at this place where the specific details of the model affect 
the speed difference v* — wjv- On the other hand, a is only 
a parameter that origin ates through the extrapolation of 
the mean-field profile (3.3) to the foremost bin region. 
We will show in the next section that the quantity a is 
a fictitious quantity, as the averag e fr ont profile deviates 
significantly from the one in Eq. ([0]) near the foremost 
bin: as we shall see, unlike the mean-field solution, it is 
not even uniformly translating. This is the reason that 
an explicit general prediction for the fronts speed beyond 
the asymptotic result obtained by Brunet and Derrida 
is hard, if not impossible, to come by. 

In passing, we note the following. It is well known from 
the analysis of uniformly translating front solutions of 
the Fisher-Kolmogorov partial differential equation that 
front solutions with v < v* are asymptotically given by 



an expression like (3^3), and that these fronts solutions 
with nodes are unstable. This does not mean, however, 
that the abov e (crude) analysis is based on an unstable 
solution ( |3.3D and therefore inconsistent. The point is 
that the expression ([0]) is only an intermediate asymp- 
totic solution, valid over some finite range of bins; just as 
in the analysis of the slow time relaxation of pulled fronts 
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in partial differential equations |17|, where such solutions 
also play a role as intermediate asymptotics, but they do 
not make the full solution unstable. 



IV. THE PROBABILISTIC DYNAMICS OF THE 
TIP: BREAKDOWN OF THE DEFINITION OF 
THE COMOVING COORDINATE £ 

We now turn to the analysis of the stochastic dynamics 
near the foremost bin, which is the region which deter- 
mines most of the front dynamics. In the light of the 
discussion of Sec. Ill B, from here onwards, we confine 



ourselves to the study of one single front realization. 

Let us assume that as the front moves in time from the 
left to the right, at some time £ = t , the bin fc is deep 
inside the saturation phase of the front. At time t > t , 
the total number of particles on the right of the ko-th bin 
is given by 



Ntot{t) 



E *fc(t) 

k>ko 



(4.1) 



For large t — to, N tot (t) grows linearly and one may define 
the asymptotic front speed vn as 



v N 



N t- 



lim 



N t ot{t) - Ntot(to) 



t 



(4.2) 



Simultaneously, the position of the foremost bin also 
shifts towards the right. For long times, the average rate 
at which the position of the foremost bin shifts towards 
the right is the same as t he f ront speed measured accord- 
ing to the definition Eq. ( |4.2| ), as otherwise, an individual 
front realization will never reach an asymptotic shape. 

Let us now examine the dynamics of the foremost bin 
in one particular realization. In Sec. III. The foremost 
bin moves towards the right by means of hops of the X 
particles. The way this diffusion takes place is as follows: 
let us imagine that in one particular realization, at a cer- 
tain time t' , the index for the foremost bin is fci, i.e., at 
time t', all the bins on the right of the fci-th bin in that 
realization are not occupied by the X particles (see Fig. 
_](a)). The diffusion of the X particles from the fci-th bin 
to the (fci + l)-th bin is not a continuous process. As a 
result, it takes some more time before the first X particle 
diffuses from the fci-th bin to the (fci + l)-th bin. Let 
us denote, by t2, the time instant at which this diffusion 
takes place (see Fig. 0(b)). Clearly, there is no exchange 
of X particles between the fci-th bin and the (fci + l)-th 
bin in the time interval £' < t < t%- During this time 
however, there can be time spans, where the number of 
the X particles in the fci-th bin may drop down to zero, 
since in the time interval t' < £ < £2, the diffusion of 
the X particles out of the fci-th bin towards its left is an 
allowed process. By definition, at time £2, the (fci + l)-th 
bin becomes the "new foremost bin" . Let us now denote, 
by t\, the time instant when the fci-th bin became the 



"new foremost bin" due to the diffusion of an X particle 
from the (fci — l)-th bin in exactly the same manner (see 
Fig. 0(c)). In this notation, therefore, £ 2 > £1, and we 
say that fci-th bin remains the foremost bin for the time 
interval A£ = £2 — t\. If we now have a series of such 
A£ values in sequence, i.e., a sequence of time values 



A£i,A£ 2 , 



, At,-, for which a bin remains the foremost 



bin, then it is easily seen that the asymptotic front speed 
is also given by 



VN 



lim j 



E **r 



(4.3) 



Put in a different way, if we denote the probability that 
a foremost bin remains the foremost bin for time A£ 
by V(At), the asymptotic front speed, according to Eq. 
.3D, is given by 



v N 



d(At) At V{At) 



(4.4) 
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FIG. 1. Snapshots of one particular realization at times t' , 
£2 and ti. The filled circles denote the X particles in different 
bins. At time £2, (fci + l)-th bin becomes the new foremost 
bin. In a similar manner, fci-th bin became the new foremost 
bin at time £1. 



Henceforth, our goal is to obtain a theoretical expres- 
sion for V(At), for given parameter values N and 7. As a 
first approach, we will make an attempt to devise a mean- 
field theory for this purpose. It is precisely at this place 







that we need to study the origin and the consequences of 
the breakdown of the definition of the comoving coordi- 
nate, £. 



A. The Stalling Phenomenon: lowest order approach 

The origin of the breakdown of the definition of the co- 
moving coordinate, £, in a mean field description is quite 
easy to understand. As can be seen from the discussion 
in the paragraph above Eq. ( |4.3| ) , the key lies in the fact 
that for the time a foremost bin remains the foremost 
bin, the front in the tip region does not move at all. We 
refer to this as the "stalling phenomenon" . During such 
stalling periods, all the dynamics is confined within the 
left of (including) the foremost bin. It is this stalling 
phenomenon that is responsible for the breakdown of the 
definition of the comoving coordinate, £ [B3|. 





the label of the foremost bin between time 
t = and t = At in an actual realization, 
e.g., in a computer simulation. 


k 


the label of the bin that attains the status of 
the foremost bin at time t = in the mean 
field theory that we describe in this section. 
Naturally, at t = 0, the density of X particles 
in it is equal to 1/N. 


km® 


the label of the bin, where the average front 
profile , extrapolated from behind, is 
equal to 1/N. 


h 


the label of the bin behind the tip, from 
which point on corrections to the profile (j)^ 
are neglected. 


k n 


The bin where becomes zero, i.e., the 
value of k where the argument of the sin func- 
tion of 4>^ becomes n. 



Table I: Summary of the various coordinate labels used in the 
paper. 

Our first step in analyzing the stalling phenomenon 
is to get back to the k and the t coordinates, but in a 
different way than we have used them so far: the fore- 
most bin, for the entire duration it remains the foremost 
bin, is indexed by an arbitrary fixed integer kt in this 
new scheme of relabelling the bin indices. The rest of 
the bins are accordingly indexed by their positions with 
respect to the fc/-th bin. Moreover, we start to count 
time (i.e., set the clock at t = 0) as soon as an X par- 
ticle diffuses into the fc/-th bin from the left and stop 
the clock just when an X particle diffuses from the kf- 
th bin to the right. This relabelling strongly resembles 
the system of comoving coordinates, hence we call it the 
"quasi-comoving coordinates" . In this formulation, the 
clock stops at time At and resets itself to zero. In this 
manner, the propagation of the front is a repetitive pro- 
cess of creating new foremost bins in intervals of At. Of 



course, it is a probabilistic process, in which the value of 
At is not fixed. 

Our mean-field theory essentially mimics the stalling 
phenomenon just as we see it in a computer simulation. 
In this theory, we also have a foremost bin, which we in- 
dex by a fixed integer k m in the quasi-comoving frame. 
In these coordinates, we describe the dynamics of the 
front by the average number of X particles in the bins. 
Between the times t = and t — At, all the dynamics 
of the front is confined to the left of (including) the k m - 
th bin. For the benefit of the reader, we summarize the 
various coordinates k used in this paper in Table I. 

The equations o f mo tion in this quasi-comoving frame, 
analogous to Eq. ( [2.5| ), in terms of the bin indices k are 
therefore given by 

^ i {N k (t))=j[{N k + 1 (t)) + (N k _ 1 (t))-2{N k (t))} 

+ (N k (t)) - i [(N 2 k (t)) - (N k (t)) } Vfc < k m , 

|W)> =7 [<JV*-i(*)> - (N k (t))} + (N k (t)) 

-J} [( N k(t)) ~ (N k (t))j iork = k m , 



and (N k ) = V k > k n 



(4.5) 



for < t < At, with the initial condition that (N km ) = 
N km — 1 at time t — 0. The angular brackets above 
denote quantities averaged over many snapshots of one 
single front realization at time t. We focus our attention 
to the region at the leading edge of the front (up to the 
k m -th bin), where the nonlinearities can be neglected so 
that the dynamics is given by 

d 

-Q t <t>k{t) =7 [<Pk + i{t) + (f>k-i(t)-2<j) k (t)} + cf> k (t), 



dt 



V k < k„ 

<t>k{t)=l [<f> k -l(t)- <f>k(t)] + M*): 



for k — k m (4.6) 



with cj) k (t) = (N k (t))/N, < t < At and 4> km = 1/N 
at time t = 0. Equation (4.5) explicitly illustrates that 
the growth of the probability ahead of the foremost bin 
is somewhat different from that behind the foremost bin 
as a result of the stalling. 

Before, we already introduced the probability T'(At) 
that the foremost bin remains the foremost one between 
the times t = and t = At. Since the foremost bin ceases 
to be the foremost one when a particle jumps out of it to 
the neighboring empty one on the right, V(t) obeys the 
differential equation 



or equivalcntly, 



l(N k Jt))V(t) 



(4.7) 



V(At)=j(N km (At)) exp 



A/ 



7 / dt(N km (t)) 



, (4. 
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satisfying the nor mali zation con dition. Clearly, as one 
can see from Eqs. (4.4) and (|4.8| ), the proper asymptotic 
speed is determined by (Nk m (t )), which in turn must 
come out of the solution of Eq. (4.6), i.e., from the effect 
of the stalling phenomenon on the leading edge of the 
front. 

The dynamics of the leading edge of the front, de- 
scribed by our mean-field theory in the previous two 
paragraphs, is a clear over-simplification. In an actual 
realization, the dynamics of the tip that governs the prob- 
ability distribution V(At) in the quasi-comoving frame, 
is quite complicated. The foremost bin has only a few 
particles, and as a consequence, the fluctuation in the 
number of particles in it plays a very significant role in 
deciding the nature of the probability distribution V(At). 
Arising out of the fluctuations, there are two noteworthy 
events that have serious consequences for the behaviour 
of V{At): (i) The creation of the new foremost bins is 
a probabilistic process, for which the time scale is char- 
acterized by l/vpf. However, if several foremost bins are 
created in a sequence relatively fast compared to the time 
scale set by 1/u/v, then one naturally expects that soon 
there would be a case when the new foremost bin would 
be created at an unusually large value of At. (ii) Ac- 
cording to our definition, in the actual realization of the 
system, the fc/-th bin remains the foremost bin between 
time t = and t = At. However, it may so happen 
that during this time, all the X particles in the fcf-th bin 
diffuse back to the left, leaving it empty for some time, 
until some other X particle hops into it, making it non- 
empty back again at a time < t < At. By the nature 
of construction, no mean-field theory can ever hope to 
capture the fullest extent of these fluctuations, and the 
one that we just presented above [that represents the ef- 
fect of the stalling phenomenon on the asymptotic speed 
selection mechanism for the front by considering V(At)], 
is no exception. Therefore, in this mean field theory that 
we described in this section, such fluctuation effects are 
completely suppressed. We will return to these fluctua- 



tion effects in Sec. IV C below, where we will make an 



attempt to estimate the effects of these fluctuations on 
T'(At). The corresponding estimates will then be used 
to improve the theoretical prediction of V{At) as well as 
to draw limits on the validity of our mean-field theory. 



B. Effect of the Stalling Phenomenon on the Front 
Shape near the Foremost Bin 



In the previous subsection, we obtained a mean field 
type expression for 'P(At) in terms of (Nk m (t)). A 
first approximation for (Nk m (t)) would be obtained from 
the solution of Eq. ( |4.5| ) above. However, in prac- 
tice the average occupation (Nk m {t)) is affected by the 
stalling effect itself. We now account for this effect in 
a self-consistent way by calculating the corrections to 
the front shape near the foremost bin. We start with 



Eqs. (4.6-4.7), and subsequently build upon the con- 
siderations of Sec. Ill, where we derived the solution 
</>(£) = Asin[zi£] exp(— z r £) at the leading edge of the 
front. 

A naive approach would be to claim that the shape 
of the lead ing edge of the front, described by the set of 
equations (4.6), is given by <ft k (t) — Asin[zi(k — vwt) + 
f3] exp[— z r (k— VNt)] for < t < At in the quasi-comoving 
frame. Notice that we have reintroduced the phase fac- 
tor f3 inside the argument of the sine function, in view 
of the fact that k can only take integral values. This 
solution of 4>k{t) would once again generate the same 
dispersion relation as in Eq. (A2). However, it is in- 
tuitively quite clear that this solution of <pk(t) cannot 
hold all the way upto k — k m , since the equations of 
motion for k < k m are different from the equation of mo- 
tion for k = k m . First of all, 4>k m {^ = 0) = V-^i which 
may not necessarily be equal to the value of the func- 
tion Asm[zi{k m — VNt) + 0\ exp[— z r (k m — v-^t)\ at time 
t = 0. Secondly, for the entire duration of < t < At 7 
the tip of the front is stationary at k m , and as a re- 
sult, the flow of particles from the left starts to accumu- 
late in the fc m -th (foremost) bin. With increasing value 
of t, bins on the left of the foremost bin get to know 
that the tip of the front has stalled, and the correla- 
tion among different bins starts to develop on the left of 
the foremost bin. As a result, an excess of particle den- 
sity beyond the corresponding "normal solution" values 
Asiri[zi(k — v^t) + (3] exp[— z r (k — v^t)} builds up on the 
left of (including) the foremost bin over time. This is 
demonstrated in Fig. ^. 



thbin 




k m th bin 
(foremost bin) 



Snapshot at time t 

FIG. 2. Snapshot of the tip of the front in a mean-field 
description at time < t < At, showing density buildup of 
X particles on and behind the foremost bin for large enough 
value of t. The dotted curve is for the "normal solution", 
Vk (t) = Asin[zi(k — fjvi) + j3\ exp[— z r (k — ujvt)]. The solid 
curve is for the actual function (f>k(t). Even though both 
(t>k°\t) an d </>/=(£) are discrete functions of k, we have drawn 
continuous curves for clarity. 



To deal with the effect of stalling phenomenon on the 
density of X particles in the bins at the tip of the front, 
which is very crucial to calculate (Nk m (t)), let us express 
4>k(t) as 



(4.9) 



where 4>f^ (t) = Asm[zi(k — vj^t) + (3] exp[— z r (k — VNt)]. 
The quantity S(j) k (t) then denotes the deviation of the 
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density of the X particles in the fc-th bin from the "nor- 
mal solution" (fyfp (t) . It takes time for the deviation 
to develop in any bin, and moreover, since such correla- 
tion effects spread diffusively, the information that the 
tip of the front has stalled at the foremost bin does not 
affect too many bins behind the foremost bin. Thus, 
it is reasonable to assume that on the left of the fore- 



most bin, there exists a bin, henceforth indexed by kb 
in this quasi-comoving coordinate (i.e., kb < k m ), where 
the magnitude of 5<f> k (t) is so small that we can impose 
t he conditio n th at 5<fi kb (t) = 0. We then substitute Eq. 
(4.9) in Eq. ( |4.6| ) and without having to worry about the 
equation of motion for 5(f>k b (t) , we obtain the equations 
of motion of the quantities 8<f) k (t) for kb < k < k m as 



— 5ct> k (t)=i [5ct> k+1 (t)-2 5(t> k (t)} + 5cf> k (t), 



for k = kb 



1. 



— 54> k (t) =7 [% + i(t) + *-i(t)-2^(t)] +64> k (t) V {k b + \)<k<k n 



h(t)} + 6Mt) - 7 [4+1 - 4 



for k = kn 



If we now denote the (k rn — fct,)-dimensional column vector [S(f> km (t), S4> km -i(t), ■ 
(4.10) becomes an inhomogeneous linear differential equation in <5<f>(£), given by 



— 6Q(t) = M6$(t) 
at 



<5$ 



p ■ 



(4.10) 



h b +i(t)} by 6$(t), then Eq. 



(4.11) 



where M is the (k m — kb) x (k m — fc(,)-dimensional tridiagonal symmetric matrix: 



M 



1-7 7 
7 I-27 

7 









7 
-27 7 







7 1 — 27 7 

7 1-27 



(4.12) 



and 5Q n 



. ^-[7(0i O 2-C+i)' '- 
straightforwardly obtained as 



, 0]. The solution of the linear inhomogeneous differential equation, Eq. (4.11), is 



<5$(t) = exp[Mi] 5<f>(t = 0) 



dt' exp[M(t - t')] 8%{t'). 



(4.13) 



To obtain the expression of (N km (t)), which is our fi- 
nal goal, we have to determine the unknowns 5Q(t = 0). 
Of these, the expression of S(j) km (t = 0) is already known 
from the fact that at time t = 0, there is exactly one X 
particle in the fc m -th bin, i.e., 



ifc™(* = 0) = ^ 



0). 



(4.14) 



The values of S(f>k(t = 0) for kb < k < k m are also quite 
easily determined when we notice that at time t = At, 
the values of 5<j) k (t = At) must reach the corresponding- 
values of 5<j)k-i{t — 0), because the average shape of the 
front repeats itself once every At time (note here that 
the repetitive character of foremost bin creation in the 
quasi-comoving frame is built in). This leads us to the 
following set of k m — kb — 1 consistency conditions 



8cj> kb+x {t = 0) 



d(At)V(At)5(f> kb+2 (At) 



5<t> km -2{t = 0) = / d{At)V(At)6<t> km -i(At) 
Jo 



(t = Q) 



d(At)V(At)S(f> km {At) 



N 



(4.15) 



The equatio n for 54> k „-i{At) is different from the other 
ones in Eq. ( [4.15 ), as it has an extra —l/N on its r.h.s. 
This is so, because the one X particle that hopped over 
to the k m -th bin at t = 0, came from the (k m — l)-th bin. 

In actuality, Eq. ( 4.15 ) should be written in terms 
of 0fe's. If we do so, then on the r.h.s. of the cor- 
responding equations, we have integrals of the form 



d(At)T(At) 



i> k (At). We have replaced these in- 



tegrals by <pu°\(t — 0). This is consistent with the fact 
that in an average sense, the underlying particle density 
field (j>^ (t) has a uniformly tran slatin g solution. The 
leftover Scj) k terms then yield Eq. ( 4.15 ). 

In terms of this formulation, the leading edge of the 
front, whose equation of motion is governed by the lin- 
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earized equation, Eq. (4.6), is divided into two parts 
|Q. In the first part, which lies on the left of (in- 
cluding) the fcf,-th bin, the solution is given by the form 
4>k(t) = Asin[zi(k m - v N t) + (3] exp[-z r (fc m - v N t)] for 
< t < At. In the second part, constituted by the bins 
indexed by k, suc h that k ^ < k < k m , the shape of front 
is given by Eqs. (4.S- Q1] ). The first part yields the lin- 
ear dispersion relation, Eq. (A2), while the second part 
yields more complicated and nonlinear relations between 
vn, z r and zi involving several other unknown quantities 
as a self-consistent set of equations. With the values of 
A, k\y and k m externally determined, if one counts the 
number of equations and the number of unknowns that 
are available at this juncture for the sele cted asym ptotic 
spee d vn , then, from Eqs. (|a|), (0), (fOl) and ([4.13 - 
4.15), it is easy to see that they involve as many un- 
knowns as the number of equations. The value of A is 
obtained by matching the mean field solution o f the bulk 
of the front, where the nonlinearities of Eq. ( f4.5| ) play 
a significant role, with the solution of the leading edge 
of the front described by the linear equations (i.e., Eq. 
(4.6)). On the other hand, obtaining the value of kt and 
k m , for a given set of parameters N and 7, is a more 
complicated process and now we address it in the next 
few paragrap hs. W e will take up these issues in further 
detail in Sec. VP as well, when we compare our theoret- 
ical results with the results obtained from the computer 
simulation. 

While it is easy to determine the foremost bin and 
hence define kf for any given realization in a computer 
simulation, the question how to obtain the values of 
k m , (3 and kb for a given set of values of N and 7, 
still remains to be answered. As a first step to an- 
swer this question, we redefine A and absorb the quan- 
tity (3 in k by a change of variable, Zik + [3 — > z;fc, 
such that in the quasi- comoving frame, (t) reduces 
to Asm[zi(k — v^t)] exp[— z r (k — v^t)]. First, this makes 
k a continuous variable as opposed to a discrete integral 
one. Secondly, the number of unknown quantities is also 
reduced from three to two, namely, to k m and fcj,. 

If we now look back at Fig. 0, and recapitulate the 
structure of the mean-field theory we presented in this 
section, we realize that the buildup of particles in the bins 
at the tip of the front due to the stalling phenomenon al- 
ways makes the curve <pk{t = 0) lie above <j>f\t = 0), 
when they are plotted against the continuous variable 
k. In our mean-field theory, <fik m {t = 0) = l/N, which 
clearly means that (f)^ (t = 0) < 1 /N and since cj)^ (t — 
0) is a monotonically decreasing function of k, this fur- 
ther implies that k m > k mo , where tf>^ (t = 0) = l/N. 

In our mean-field theory, what is the numerical value 
of (k m — k mQ ), the distance between the bin, where the 
lowest order approximation cj)^ reaches the values l/N 
and the bin, where the actual average profile <fi reaches 
this value? For arbitrary values of N and 7, this is not 
an easy question to answer. 



To check our theory, in this paper we confine ourselves 
mostly to the case of 7 = growth rate= 1, as it is the 
most illustrative case to demonstrates the multiple facets 
of fluctuating front propagation. For a part of the analy- 
sis, we also consider the 7 = 0.1 case. For such values of 
7, i.e., if 7 is too small (7 <§; 1), or not too large (7 ~ 1), 
the only information that we have at our disposal to ob- 
tain the value of the continuous parameter k m , is the 
fact that k m > k mo . For such values of 7, therefore, the 
only remaining way to generate the V(At) curve is to use 
trial values of k m , for k rn > k rno in an iterative manner 
jyj [recall that the value of k m is needed for the initial 
condition, Eq. (4.14)]. For such values of 7, the use of 
the trial values of k m to generate V{At) also requires the 
value of k m — k^ as an external parameter, which can be 
chosen to be a few, say ~ 4 (of course, this number can 
be increased to obtain higher degree of accuracy for the 
6cf>k(t = 0) values). We will take up further details about 
it in Sec. V. However, before that, we next discuss two 
additional fluctuation effects that have important conse- 
quences on the V{At) curve. We also mention here that 
we have explored the possibility of a relation between 
kf, obtained from computer simulation results, and k m , 
but due to the fact that kt has stochastic fluctuations in 
time, such a relation does not exist. 



C. Additional Fluctuation Effects 

Having described the mean-field theory, we are now in 
a position to assess its accuracy or validity for the prob- 
ability distribution V(At) that it generates, before we 
start to look for numerical confirmation. At the end of 



Sec. [VA, we have mentioned that the fluctuation of the 



number of X particles in the foremost bin plays a very 
significant role in deciding the nature of V{At). Such 
fluctuations are not captured in our mean field theory, 
which simply assumes that the number of X particles in 
the foremost bin at t = is 1 and afterwards the number 
of the X particles in it increases through the proc ess of 
a mean growth. In particular, at the end of Sec. IV A 
we have described two kinds of events that, we now ar- 
gue, affect the nature of V(At) for large values of At, 
compared to the time scale set by 1/ujv- 



1. Few Foremost Bins Are Created Too Fast in a Sequence 

The first of these events is that if a few of the new fore- 
most bins are created relatively fast in a row, then soon 
there would be a case of a new foremost bin creation 
that takes an unusually long time. Naturally, this gives 
V{At) a higher value than what our mean-field theory 
does for large values of At. The reason for this is quite 
simple: the mean growth of the number of X particles 
in the foremost bin is exponential in time, which would 
indicate that if one describes the growth of the number of 
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X particles in the foremost bin simply by mean growth, 
then the probability distribution V(At) decreases very 
rapidly for large At, and clearly that fails to describe the 
slow decay of V{At) for large At arising out of this event. 

Unfortunately there is no way to estimate the effect 
of this event within the scope of any mean-field theory, 
since by its sheer nature, it can only be described by the 
multi-time correlation functions of the times required for 
sequential creations of new foremost bins. For this rea- 
son, we call this event "correlated diffusion event" for 
later reference. But the physical effect of it can be ex- 
pressed in a slightly different manner which is more con- 
ducive for understanding the conditions of applicability 
of our mean-field theory. In our mean-field theoretical 
description, before a new foremost bin is created, the 
shape of the fro nt is always the same mean shape, de- 
scribed by Eqs. (ff.5|). On the other hand, if a few of the 
new foremost bins are created relatively fast in a row, 
the leading edge of the front gets more and more elon- 
gated while the number of particles inside the bins in 
the leading edge does not get a chance to grow accord- 
ingly. Thus, this event creates significant deviation for 
the actual front shape from the front shape described by 
our mean-field theory. The magnitude of this deviation, 
measured by subtracting the mean-field density of the X 
particles from the actual density of X particles inside the 
bins at the leading edge of the front, is obviously nega- 
tive. If we combine this argument with the fact that on 
an average, the probability of a new foremost bin cre- 
ation increases with the increasing number of X particles 
in the foremost bin, then it is easy to realize that after a 
sequence of such fast creations of new foremost bins, the 
front needs to replenish the number of X particles in the 
leading edge before another new foremost bin is created. 
It is this replenishing process which is responsible for the 
next new foremost bin creation at a relatively long time. 

It is now intuitively clear that in terms of the front 
shape, the larger the deviation such an event causes, the 
more V(At) will be affected for large values of At. Based 
on this, we now argue that for a fixed value of N, such 
an event does not affect the large At behaviour of V(At) 
curve for large values of 7 as much as it does for small 
values of 7. To reach this conclusion, one simply needs 
to observe the following: the mean shape and the cor- 
responding density of the X particles in the bins at the 
leading edge of the front is characterized by z r and Zi, 
and for small 7, the values of z r and Zj is large and vice 



versa (as ~ z r ~ 7 ' , see Eq. (AE)). For large 7, 
therefore, for the mean shape of the front, the leading 
edge is already quite elongated and the density of the X 
particles at the tip of the front is quite small, compared to 
their small 7 values. As a result, for small 7, the magni- 
tude of the deviation from the mean front shape, caused 
by such an event is much larger, and since the growth 
rate is always unity for all 7, it takes a much longer time 
to replenish the density of X particles for small 7 than 
for large 7 values. 



2. Particles in the Foremost Bin Jumps Back 



The other kind of fluctuation effect has to do with the 
fact that albeit according to our definition, the foremost 
bin remains the foremost one until time t — At, it may 
so happen that at some nonzero value of t, all the X 
particles diffuse back to the left leaving the foremost bin 
empty for some time and then another X particle diffuses 
into the foremost bin from the left, making it non-empty 
again at a finite value of t, say at t = to. Clearly, this 
event is much more unlikely to take place once the num- 
ber of X particles in the foremost bin has grown, since in 
that case, all the X particles in the foremost bin have to 
diffuse back to the left. Essentially, this event is there- 
fore restricted to the following sequence: (a) starting 
t = 0, the foremost bin remains occupied by a single X 
particle for sometime, (b) this X particle then diffuses 
back to the left leaving the foremost bin empty, until (c) 
another X particle diffuses into the foremost bin, making 
it non-empty again at t = t$. Of course, the value of to 
is not fixed and it is chosen probabilistically. For later 
reference, we call this "the vacant foremost bin event" , 
and this event is much more likely to take place for large 
diffusion coefficient 7. 

Based on the picture described in (a)-(c) above in the 
previous paragraph, we can now make a quantitative es- 
timate of this particular event and accordingly correct 
the expression of V(At). One simply has to realize that 
if this event takes place, then the time at which the the- 
ory for mean growth of X particles in the foremost bin 
(with exactly one X particle to start with) can be ap- 
plied in this quasi-comoving frame, shifts from t = to 
t = to- However, we also need to obtain an estimate for 
the value of to. This can be obtained using the following 
argument: if in step (b), the only X particle in the fore- 
most bin had diffused to the right, instead of diffusing to 
the left, it would have been a case of a new foremost bin 
creation, the time scale for which is set by 1/u/v- Since 
the probability of this single X particle in the foremost 
bin to jump to the right is the same as the probability 
of it to jump to the left, we can also say that the time 
it takes for the X particle in the foremost bin to diffuse 
back to the left takes approximately a time 1/v/v start- 
ing t = 0. Similarly, step (c) is exactly the same step 
as a "new foremost bin creation". Hence, after step (b) 
is over, it takes a further 1/wjy time [M for another X 
particle to diffuse from the left into the foremost bin. To- 
gether, these two events make to « 2/vn W^ i an d this 
argument illustrates that this event affects the behaviour 
of V(At) only for At > 2/vn- Having neglected the effect 
of the correlated diffusion events on the 'P(At) curve for 
t > 2/vn (for which we have no theoretical estimate any- 
way), if we now claim that for all At values greater than 
2/vn, the population of the X particles in the foremost 
bin is described by Eqs. ( |4.5| ), but with the condition 
that (N km (t = t )) = N km (t = t ) = l/N, as opposed to 
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having (N km (t = 0)) = N km (f = 0) = 1/N, then we can 
still incorporate the effect of this event (that arise out 
of fluctuations) within the scope of the mean-field theory 
that we described in this section. If this procedure is cor- 
rect, then while comparing the theoretical V(At) curve 
with the V(At) curve obtained from the simulations, one 
would notice that for large values of At, this procedure 
underestimates the magnitude of V(At). Hitherto, this 
underestimation then would be an indication of the effect 
of the correlated diffusion events on V{At). We will re- 
turn to these points once again in the next section, where 
we seek numerical confirmation of our theory pSj . 



D. Summary of the Status of the Present Approach 
and Additional Observations 



The discussion above completes the theoretical formu- 
lation for the asymptotic speed selection of the front. 
Before we discuss how given values of N and 7 would 
generate the corresponding values of vn from our mean- 
field theory described in this section, we summerize our 
claims here and make a number of additional observa- 
tions. 

1 ) Based upon the microscopic description of the front 
movement, we have formulated a mean-field theory that 
describes, on a lattice, the front propagation as a se- 
quence of "halt-and-go" process. In this way of looking 
at the front propagation, essentially the number of X par- 
ticles at the tip of the front determines the asymptotic 
speed of the front. Since the number of X particles at 
the tip of the front are rather few, the fluctuations in 
the number of X particles at the tip of the front affects 
the asymptotic speed of the front in a strong manner. 
Part of the fluctuation effects can be estimated within 
the scope of this mean- field theory itself. The other part, 
for which the fluctuations can only be studied by means 
of a multi-time correlation functions, is expected to af- 
fect the accuracy of our theory much more for small 7 
than for large 7 values. Therefore, overall, in terms of 
numerical confirmation, one can expect to find a greater 
accuracy for large values of 7. 

Moreover, for large 7, the discreteness of the lattice 
effects are suppressed, and therefore, for a given value 
of N , one would expect that the relative correction for 
the asymptotic front speed, (v* — vm)/v* , must become 
small. 

2) There are two important aspects that one must take 
notice of. First, in a mean field description that incor- 
porates the effect of the stalling phenomenon, we have 
demonstrated from the microscopic dynamics that there 
exists a cutoff of particle density , which is expressed by 
the fact that in this mean field description, 4>k{t) = 
for k > k m against a finite value at k = k m . Secondly, 
we have also demonstrated that the quantity a in Sec. 
Ill is indeed an effective quantity only, as the soluti on o f 
the linearized equation of the front, given by Eq. (3.3), 



is not valid near the foremost bin, and the fact that 
the asymptotic speed selection mechanism arises from a 
proper probabilistic description of the tip of the front. 

3) For very small values of 7, we have previously no- 
ticed that the correlated diffusion event plays a very dom- 
inant role that no mean-field theory can ever generate, 
so we should leave the 7 <§; 1 case outside the purview 
of our mean-field theory (we will demonstrate this in the 
next section). 

4) To judge the appropriateness of our mean-field the- 
ory, as far as the generation of the numerical value of the 
asymptotic front speed for given values of N and 7 is con- 
cerned, we make the following observations: (i) the case 
of 7 <C 1 cannot be studied in terms of a mean-field the- 
ory, (ii) the case of 7 ~ 1 needs a trial value of k m > k„ lQ 
and the use of a recursive feedback mechanism to gen- 
erate the V(Af) curve, and (Hi) we still need the values 
of A, z r and Zi, which can be obtained only from the 
simulation data for given values of N and 7. In view of 
these points, it is clear that this theory is unable to make 
a definitive predictions for vn, without any assistance 
from the compute r simulations whatsoever. Moreover, 
the Eqs. (4.4- 4.15 ), which one needs to solve to generate 
the V(At) curve, are highly nonlinear equations, hence, 
this theory can only hope to show consistency with the 
results of the computer simulations, as opposed to pro- 
duce a numerical value of vn which is then subsequently 
confirmed by the computer simulations. 

5) Finally, we note that unlike Eq. ( |3.4[ ), this theory 
does not make the effect of the value of N on the asymp- 
totic front speed explicit. However, it is natural to expect 
that the effect of stalling of the front and the associated 
particle density buildup at the tip of the front on the 
front shape and speed would become less and less for in- 
creasing N. This would reflect in the comparison of the 

&4>k{t = 0) values against the corresponding <p£\t = 0) 
values in the bins at the very tip of the front. We would 



return to this point in Sec. VD2 



V. TEST OF THE THEORY AGAINST 
COMPUTER SIMULATIONS 

We now check our theory, as it has been presented in 
Sees. Ill and IV, against the results of the computer 
simulations. There have been quite a few aspects of the 
theory that we have presented in Sees. Ill and IV; and for 
a given set of values of N and 7, testing all these aspects 
of our theory is not a short and easy process. To explain 
how we do the simulations, obtain vm and A, and check 
the front shape, we choose one particular set of N and 7 
values, namely N — 10 4 and 7=1. We then use these 
methods to obtain the simulation data for three other 
values of N, namely N = 10 2 , N = 10 3 and N = 10 5 , 
keeping the value of 7 fixed at 1. Based on this scheme, 
this section is divided into five subsections. In Sec. VA, 



we present the simulation algorithm and obtain vn for 



15 



7 = 1 and N — 10 . In Sec. VB, we summarize the 
method to calculate A, and subsequently obtain its va lue 
for 7 = 1 and N — 10 4 using the results of Sec. VA. In 
the Sec. V C , we contrast the simulation results of Sees . 
VA| and |VB| with the theory of Sec. III. In Sec. |VD|, 



we test our theoretical predictions for V(At) against the 
computer simulation results for N = 10 4 10 2 , 10 3 and 
10 5 , and 7 = 1 (in that order). Moreover, in Sec. IV, 
we have conjectured that the mean-field theory mimick- 
ing the stalling phenomenon would be less success ful for 
small values of 7. We verify this conjecture in Sec. VD 2 



by means of a relative comparison of the theoretical and 
simulation V(At) curves for 7 = 0.1 and N = 10 4 . We 
also remind the reader that in Sees. VA-VC, fc and t 



respectively denote the laboratory bin co-ordinate and 
actual physical time (and therefore they do not relate to 
the quasi-comoving co-ordinates or the resetting of clocks 
that requires < t < At). 



choose another set of two random numbers, n and r2, in 
[0, 1). From the numerical value of r-y and the fact that 
the probability of a transition taking place in the fc-th 
bin is given by Wk/Wtot, we determine the index of the 
bin where the transition takes place. Similarly, the par- 
ticular transition in the fc-th bin is determined from the 
numerical value of T2 and considering the probabilities of 
different kinds of transitions in the fc-th bin: 



probability of a diffusion to the right = 
probability of a diffusion to the left = 



probability of breeding a new X particle = , and 



probability of annihilating an X particle = 



N k (N k -l) 



NW k 



(5.5) 



A. Computer Simulation Algorithm 



Our algorithm for carrying out the computer simula- 
tions is the same as it has been described in Ref . || . The 
starting density profile of the X particles is a step func- 
tion, given by <j)k{t = 0) = [ 1 — 0(fc — fco) ], for some fco. 
The simulation algorithm consists of a repetitive itera- 
tions of two basic steps: 

(i) Let us assume that at any time i, the configuration 
of the system is given by (N±, N2,..., N k >), for some k' . 
The total rate of possible transitions, W k , for the N k 
number of X particles in the k-th bin are the sum of 2N k 
diffusions, creation of N k new X particles and annihila- 
tion of N k (N k — 1)/AT number of X particles, i.e., 



W k = 2jN k + N k + 



N k (N k 



1 



N 



(5.1) 



The total rate of transition, Wtot, for all the bins is there- 
fore 



Wt, 



k' 

E 

k=l 



(5.2) 



Starting at time t, the probability of no transition hap- 
pening for an interval r is given by 



p(r) = exp(- Wtot t) . 



(5.3) 



Before any transition takes place, a random number ro 
is chosen within the interval [0, 1). The time r that one 
needs to wait before any transition happens is then de- 
termined as 



1 



lnr 



(5.4) 



(ii) With the time r for a transition at our disposal, the 
bin where the transition takes place and the specific tran- 
sition in that bin must also be determined. To do so, we 
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FIG. 3. To illustrate that the front reaches its steady 
state shape before t = 200, the plot of 4>k{t) vs - t 
for five different values of t spaced at regular intervals, 
t = 200, t = 250, t = 300, t = 350 and t = 400, are shown 
above. 



Once the transition is determined, the configuration of 
the system is subsequently updated. However, any X 
particle diffusing from the first bin (i.e., k = 1) towards 
the left is immediately replaced. 

In this subsection, we focus on one particular set of val- 
ues of 7 and N, namely, 7=1 and N = 10 4 . The value 
of fco for the initial density profile of the X particles is 
chosen to be 50. Starting at t = 0, we let this initial 
profile evolve in time. To obtain the random numbers, 
we use the random number generator drand48 provided in 
the standard C library functions with the initial seed |42] 
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s = 123456. It turns out that to a very good approxima- 
tion, the front shape reaches a steady state somewhere 
before t = 200. The front shapes from t = 200 to t = 400 
is shown in Fig. [| as an illustration. For measurement 
of the asymptotic quantities, therefore, we take t = 200 
as our starting point. 



2.75 



2.25 
VP 



1.75 



1.25 




1.964 



B. The coefficient A as a reflection of the nonlinear 
front behavior 



The quantity A has been introduced to solv e for the 
linear difference-differential equation, Eq. (3.1). Its nu- 
merical value, however, cannot be determined from the 
linear equation, since any value of A satisfies it. To de- 
termine the value of A, therefore, one needs to solve the 
full nonlinear difference-differential equation, Eq. (|2.6|) , 
expressed in terms of the comoving coordinate, £. This 
is done, together with the associated values of the real 
and imaginary part of z, in Appendix |bJ for 7=1 and 
AT = 10 4 , we find 



200 300 400 500 600 700 



A = 0.961 ± 0.012 . 



C. Numerical Test of the Predictions for 
Consistency of Front Shape and Speed 



(5.9) 



FIG. 4. Values of v N (tj) for 200 < tj < 700 and 
j = 1,2, ...,980. As expected, the fluctuations in iijv(tj) 
die out for large values of j. 



To calculate the asymptotic speed of the front, we mea- 
sure the At values for creating new foremost bins after 
t = 200 till t = 700. We find that altogether there 
are 980 different At values in this time interval. As- 
suming that the j-th value of At takes place at time tj 
(j = 1, 2, . . . , 980), we define the j-th cumulative average 
of the At values as 



(At), 



1 1 

7 E At r > 



(5.6) 



j'=l 



which subsequently allows us to define the speed at time 



tj as 



(At), 



(5.7) 



Naturally, for small values of j, the values of VN(tj) fluc- 
tuate, but as j becomes large, the fluctuations die out 
and vjs[{tj) approaches vn- The plot of VN(tj) vs. tj 
is shown in Fig. | for j = 1,2,..., 980, t x = 200.562 
and tgso = 699.271. We notice from the plot that the 
fluctuations in Wjv(tj') are really small for tj > 500. The 
VN(tj) values for tj > 500, therefore, allows us t o set the 
error bar on the measurement of in Eq. (5.S), and we 
obtain 



VN 



1.964 ± 0.006. 



(5.8) 



Equipped with the value of A and the simulation data, 
we are now in a position to contrast the result of Sec. Ill 
with the simulation results. The purpose of this subsec- 
tion is twofold: first, we demonstrate that the theoretical 
shape of the front generated in Sec. VB for 7 = 1 and 
N = 10 4 agrees very well with the shape of the front ob- 
tained from the simulation data by taking snapshots at 
arbitrary times. Secondly, we demonstrate that there are 
significant differenc es i n the two values of v* — Uj y, o ne 
obtained from Eq. (pTq), and the other from Eq. (pll). 



1.5 



• Front shape at t= 200 
□ Front shape at t= 280 
o Front shape at t = 360 
Front shape obtained from theory 




452 



FIG. 5. Theoretical shape of the front represented by the 
solid line and the shape of the front obtained from computer 
simulations at three discrete times, t — 200, t — 280 and 
t = 360, represented by three different symbols. 
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To compare the shapes of the front obtained from the 
theory and the simulations, we choose to take snapshots 
at three discrete times, t — 200, t = 280 and t = 360. 
Below we present the plot that contains the front shapes 
for these times obtained from the simulations and the 
theory. The comoving coordinate £ for the bins is chosen 
in a way such that it coincides with the laboratory frame 
position of the front at t — 200. As we can see from the 
graph below, the collapse of the data for t = 200, t = 280 
and t = 360 is very good and the solid line represent- 
ing the theoretical prediction is almost indistinguishable 
from the simulation data. The theoretical curve is first 
generated using the method described in Sec. VB with 
«w(i98o) as the asymptotic front speed, the correspond- 
ing values of Zi, z r , Zi and A of Eqs. (B6) and (fT^), and 



then having it shifted to coincide with the laboratory 
frame position of the front at t = 200. 

We now return to the re sult of Sec. Ill, and denote 
v* — vn obtained from Eq. (3J3) by Av asvrn p. For 7=1, 
Eq. (|^) yields z = 0.9071032.. and v* = 2.07344. Us- 
ing these values for N = 10 4 , we obtain 



Av a 



>///// j> 



= 0.152024. 



(5.10) 



On the other hand, using i>jv(i98o) f° r v n, the value 
Av s im = v* — vn comes out to be 



Vjv(^98o) 



0.11 



(5.11) 



which implies that the asymptotic estimate Av asyrn p is 
about 38% larger than Av sym from the computer simu- 
lations. 

These results clearly indicate that for large but not ex- 
tremely large values of N (N = 10 4 here), there is much 
more to the story than v* — ujv being simply oc In -2 N. 
The theory presented in Sec. Ill does capture the essen- 
tials, and it would have been good enough to generate 
appropriate numbers for v* — u/v, if one could obtain the 
value of a externally. However, in view o f th e fact that 
the uniformly translating solution of Eq. (2.6) cannot be 
extended all the way up to the foremost bin, the quantity 
a is a fictitious and simply an effective quantity (already 
mentioned in Sec. III). Therefore, it is not possible to 
obtain the numerical value of a from computer simula- 
tion results or from any theoretical estimate. Besides, the 
theory of Sec. Ill completely overlooks the microscopic 
intricacies at the tip of the front, and hence, it should be 
regarded as an effective theory. 



D. Numerical Test of the Theory for V(At) 

In this subsection, we seek the numerical test of our 
theory presented in Sec. IV. We carry out this task in 
two steps. In the first step, we check most of the as- 
pects of the theory for N — 10 4 and 7 = 1, where we 
describe the method for obtaining V(At). Subsequently, 
in the second step, we check the predictions of our the- 
ory for N — 10 2 , 10 3 and 10 5 , keeping the value of 7 



fixed at 1. Notice that comparing probability distribu- 
tions allows us to verify more detailed representations of 
the actual forward movements of the foremost bin against 
comparing only the asymptotic front speed vn, w hich is 
the inverse of the first moment ofV(At) [see Eq. ^2 )J. 
We should note that in view of the strong nonlinearity of 
the self-consistent theory of Sec. IV, we will have to use 
the values of A, z r , vn and z, obtained from computer 
simulations to generate the V(At) curve, and then ob- 
tain the theoretical value of vn- This process therefore 
becomes a self-consistency check of our theory of Sec. 
IV, as opposed to a verification of its predictions. More- 
over, we do not compare the V{At) curves directly. This 
is for a very simple reason: namely that the expression 
for V{At) in Eq. (|J) involves (N km (At)) as a coeffi- 
cient. In an actual computer simulation, this quantity 
fluctuates wildly, and hence, generating a histogram to 
obtain the probability distribution V(At) from computer 
simulations proves to be difficult. Instead, we compare 
the "cumulative probability distribution" curves P(At), 
which is defined as the probability of a new foremost bin 
creation happening at time t > At. From a theoretical 
point of view, the expression of P(At) can be found easily 
from Eq. (O) as 



poo 

P(At) = / dt'V{t') 

J At 



I At 

exp 



pAt 

- 7 / dt(N km (t)) 
Jo 



(5.12) 



It turns out the P(At) histogram generated from the 
computer simulation results is not noticeably affected by 
fluctuations, which makes its comparison with the P(At) 
curve generated from our theory much simpler. 



1. The Case 0/7 = 1 and N = 10 4 

The P(At) curve from the computer simulations are 
generated in the following way: by definition P(0) = 1. 
For N = 10 4 and 7=1, there are 980 values of At. First, 
these are arranged in an increasing order of magnitude, 
Ati , At2 , • ■ ■ , Atgso , and then in the corresponding values 
of P(At) are obtained as 



P(At j+1 ) = P(Atj) - 



1 

980 



(5.13) 



for j = 2,3, ...,980 with the initial condition that 
P(Aii) = 1 - 1/980. The corresponding P(At) vs. At 
plot is shown in Fig. ^| by open circles. 

To generate the corresponding theoretical cumulative 
probability distribution, we proceed in the following way. 

In a coordinate system, where the function 0^ (t = 0) = 
at k — 7r / Zi [which allows us to use the values of A, z r , 
vn and z% of Eq. ( |B6| ) j, we wo rk out the whole machin- 
ery described by Eqs. (4.5- 1.15| ) neglecting the fluctuation 
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effects described in Sec. IV C. This process requires the 
value of k m — kb, i.e., the number of bins at the tip of 
the front where the buildup of particle density is signifi- 



cant, as well as the value of fc„ 



as external inputs, 



and we choose k„ 



kb = 4 for this purpose |40| (we re- 

The 



fer back to Table I for the definitions of k mo etc 
calculation of the value of k m — k mo and the generation 
of the P(At) curve are carried out self-consistently and 
hence by iteration, using the recursive feedback method 
However, to generate the P{At) curve for any guess 
value of k m — k mo , one still needs to have the values of 
6&(t — 0) for the k m — k^ bins at the tip as a starting 
point [see Eq. ( 4.13) )] . At the same time, we notice that 
6&(t = 0) can only be determined once the probability 
distribution P(At) is obtained. We choose to address 
this problem the following wa y: fo r any guess value of 



we start with Eq. ( ]f.l4|) and 5cj> k (t = 0) = 



for the rest of the k m — kb bins. Keeping k m — k mo fixed, 
we then generate the corresponding V(At) curve and ob- 
tain the S(f>k(t = 0) v alues for the rest of the k m — kb 
bins using Eq. ( [4.15 ). We keep repeating this process 
until we converge in terms of the S(f> k (t — 0) values, i.e., 
when the recursive correction to the values of 5<j) k (t = 0) 
becomes less than 10% of the Scj) k (t = 0) values at the 
previous step in the recursion. Once this point is reached 
for a value of k m — k mo , we then compare the theoretical 
cumulative probability distribution P(At) with Fig. || 
above to decide upon the next guess value of k m — k mo 
in the recursive feedback method. 



Computer simulation data 

Fully consistent theory of Sec. IV 
With 5<|> t a=0) = for k * k m 
in the theory of Sec. IV 




2.5 



FIG. 6. The cumulative probability distribution P(At) as 
a function of At for N = 10 4 and 7 = 1. 



For N = 10 4 and 7 = 1, the value of k„ 



- k mo turns out 



to be k„ 



1.1431. We present the corresponding 



theoretical cumulative probability curves in Fig. |6|. The 
solid line in Fig. re presents the fully consistent solu- 
tion of Eqs. (p~5|-4.15), while the dashed line represents 
the theoretical cumulative probability curve obtained by 



solving Eqs. (pj-|4.15D with 8<j) k (t = 0) = for k ^ k m . 
The fact that the fully consistent solution curve matches 
the computer simulation one much better than the naive 
approximation where all the S(j) k corrections behind the 
foremost bin are ignored is a strong indication of how 
significantly the buildup of particle densities in the bins 
behind the foremost one affects the property of P(At). 



Computer simulation data 
Improved theoretical curve 




At 

FIG. 7. The theoretical curve, which includes the effect of 
the vacant foremost bin events, and the simulation data for 
the cumulative probability distribution P(At) are presented 
above for N = 10 4 and 7 = 1. 



An examination of Fig. |6| immediately reveals that 
the agreement between the P(At) curve generated by 
the fully consistent theory of Sec. IV and the computer 
simulation is extremely accurate up to about At = 0.9. 
However, the theory is unable to capture the "tail" of 
the P(At) curve for large At. Analysis of our data 
shows th at this is due to the fluctuation effects discussed 
in Sec. [VC. As mentioned there, correlated diffu- 
sion events are not captured in this theor y. However, 



we can follow the argument of Sec. [V C 2 to take into 
account the effect of the vacant foremost bin events on 
the P(At) curve for large At values: we assume that all 
cases of At > to w 2/vn are due to the vacant foremost 
bin events. This means that for At < to, the P(At) 
curve is given by the solid line in Fig. ^, but from to 
onwards, the P{At) curve must be generated from the 
mean-field dynamics of the tip described in Sec. IV, 
with the same value of k m , but with the initial condi- 
tion (N km (t = t )) = N km (t = t ) = 1/N, as opposed 
to (N km (t = 0)) = N km (t = 0) = 1/N. With the value 
of k m already determined in this subsection, the corre- 
sponding equations, Eqs. (p~5-4.15), are easy to solve 
self-consistently as before. From this analysis, we obtain 
the behaviour of P(At) for At > to, having noticed that 
P(At — > to+) must be the same as the value obtained 
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from the solid line in Fig. ^ at of At = to, i.e., 0.06034. 
We present the final theoretical P (At) curve together 
with the simulation data in Fig. [7]. Notice that this pro- 
cess introduces a finite discontinuity in the density of the 
X particles inside the foremost bin at At = to, since for 
At < to, the density of the X particles in the foremost 
bin is obtained from a fully consistent theory of Sec. IV, 
while at At = to, it is set equal to 1/iV manually. As a 
consequence, the theoretical P(At) curve in Fig. [?] has a 
slope discontinuity at At = to. 

In Fig. ^, the improved theoretical curve follows the 
curve reasonably well, but it still lies below the simula- 
tion data points for At > 2/vn, as it should be. This 
discrepancy gives us an idea about the effect of the cor- 
related diffusion events o n th e P(At) curve that we could 
not estimate. Using Eq. (5.7) to calculate the front speed 
from the theoretical curve in Fig. 0, we obtain 

v N (theoretical) = 1.98882. (5.14) 



This is about 0.024 higher than the asymptotic front 
speed measured by the computer simulation [see Eq. 
(|5.8D ], in agreement with the fact that the theoretical 
curve for -P(At) slightly underestimates the simulation 
one for At > 2/vjy. 



2. The Cases of N = 10 2 , 10 3 and 10 5 , with 7 = 1 



We now further test our theory for N = 10 2 , 10 3 and 
10 5 , keeping the value of 7 fixed at 1. The values of 
vn (simulation), z r , Zi and A in Table II below. The cor- 
responding P(At) vs. At graphs, which are the analogs 
of the graph in Fig. fj], have been plotted together in 
Fig. ||. Table II presents the theoretical values of vn 
that are calculated using t hese P(At) vs. At graphs, 
and predicted vn from Eq. (|3.6|). 



N 


vn (simulation) 




Zi 


A 


v n (theoretical) 


v N [Eq. (3.6)] 


10 2 


1.778 


0.8217 


0.436 


0.8836 


1.808 


1.465 


10 3 


1.901 


0.8586 


0.3313 


0.9042 


1.899 


1.803 


10 5 


2.001 


0.8885 


0.2654 


1.0714 


2.057 


1.976 



Table II: The vn (simulation) , z r , Zi, A, t)jv (theoretical) and tijv[Eq. 
values are calculated from the theoretical curves of Fig. ^| 



Computer simulation data 

Fully consistent theory of Sec. IV 



& 
I — 



P(At) 




values for 7 = 1, and N = 10 2 , 10 3 and 10 5 . The tijv (theoretical) 



FIG. 8. The combined theoretical curve and the simula- 
tion data for the cumulative probability distributions P(At) 
vs. At for N = 10 2 , 10 3 and 10 5 , and 7 = 1. The curves for 
the latter two are shifted upwards for clarity. 

Notice that as N decreases, according to Table II, the 
value of z r decreases more and more from its N — ► 00 
limit zq, while zi increases. This is an illustration of how 
the non-mean field effects become increasingly important 
behind the tip region. We should also note two more 
points about Fig. |[ (i) in the absence of any estimate of 
the correlated diffusion events for At > 2/vjy, the theo- 



retical curves lie below the simulation data (although for 
N = 10 2 and 10 3 , it is not so clearly discernible), and (ii) 
the agreement between the theoretical P (At) curve and 
the simulation one for N = 10 5 may appear to be worse 
than the corresponding ones for N — 10 2 , 10 3 and 10 4 , 
but this may due to the fact that we have had to con- 
tinuously cut off particles from the saturation region of 
the front on the left to obtain the stochastic simulation 
data for N — 10 5 within a reasonable computer time. 
We have found that the shape of the -P(At) histogram 
obtained from the simulation gets slightly modified de- 
pending on how this subtraction in carried out, specially 
in the At > 2/vn region. 

We now return to (the issue raised in point 5) of Sec. 



IV D) the increased importance of the stalling effects and 

00 asymptotic theory for 
shows the comparison 
0) = 



the deviations from the N 
decreasing values of N . Figure |_ 
between the actual particle numbers, (SNk)(t 

N5<j> k (t = 0) and the (N[ 0) )(t = 0) = V^ 0) (t = 0) val- 
ues for four foremost bins, i.e., for k — k m , fc TO _i, k m — 2 
and fc m _3 (note that for the sake of clarity, we have omit- 
ted the angular brackets for notations in Fig. ||). These 
values have been obtained self-consistently, while gener- 
ating the theoretical curves of Figs. and ||. As ex- 
pected, it is clear that the 6Nk(t = 0) values are playing 
less and less role for increasing value of N. There are cou- 
ple of more points that one must take notice of. First, as 
can be seen from Table II, A increases with N, but only 
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by a small amount. Secondly, it is also clear from Fig. 
g that w.r.t. the fc n -th bin [where (N^) vanishes], the 
position of the fe ro -th bin (where (N^ + SN)(t = 0) = 1) 
shifts gradually towards the left for increasing TV (see 
Table I for the definition of k n ). All these together elu- 
cidate that for not very large values of TV, the actual TV 
dependence of the front speed vn is a much more com- 
plicated story than simply the 1 / In 2 TV relaxation to v* 
oi vn- From the trend of the gradual left-shifting of k n 



(w.r.t. k n ) and the gradual unimportance of the role of 
8Nk(t = 0) values compared to their N^(t = 0), it is con- 
ceivable that for extremely large values of TV, k m — > k„ lQ 
and 5Nk(t = 0) — * 0, and it is this limit where the cutoff 
(at = 1/TV) picture in Ref. |l becomes applicable. 
In this sense, the theory of Sec. IV is complementary to 
that of Ref. |)| , as together they span the whole range of 
TV values, from reasonably small to very large. 
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10 4 and 10 5 



9. Comparison between the {SN^ft = 0) and the {N^)(t — 0) values for four foremost bins and for TV = 10 2 , 10 3 



The angular brackets for the notations have been omitted in the figure for clarity. Note that as N increases, the 

r(0) 



corrections (SNk)(t = 0) compared to the (N^ '){t = 0) profile become less and less important. 



r 



E. The Case of 7 = 0.1 and TV = 10 4 

We now investigate the claim made in Sec. IV that 
the correlated diffusion events affect the 'P(At) curve so 
severely for low diffusion coefficients that our approach 
fails badly, by comparing the theoretical P(At) curve 
with the simulation one, for 7 = 0.1 and TV = 10 . 

We present the two curves in Fig. |l^. The asymptotic 
speed for the corresponding pulled front, v*, for this set 
of parameter values is given by v* = 0.7754 and the sim- 
ulation results yield = 0.698. 

The theoretical curve of Fig. is analogous to that 



of Fig. U represented by the solid line, and it is obtained 
by means of a fully consistent theory of Sec. IV. Notice 
that the agreement between the theory and the computer 
simulation results is not good beyond At « 1/vat. It is 
also obvious that an attempt to incorporate the effect of 
the vacant foremost bin events does not do any improve- 
ments in this case, since the value of -P(Ai) is almost 
zero for At > 2/vn- This is very much expected and a 
careful examination of the simulation data also reveals 
that the vacant foremost events do not occur at all dur- 
ing the front speed measurement times between t = 200 
and t = 700 (see Eq. (5.6) and the paragraph above it). 
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Altogether, this fits very nicely in the consistent picture 
that we have put forward so far, which simply indicates 
that the entire discrepancy between the theory and the 
computer simulation in Fig. |To| is solely due to the cor- 
related diffusion events. 



1.2 



Computer simulation data 

Fully consistent theory of Sec. IV 




FIG. 10. The cumulative probability distribution P(At) as 
a function of At for N = 10 4 and 7 = 0.1. 



VI. THE CASE OF A FINITE NUMBER OF Y 
PARTICLES ON THE LATTICE SITES 



each lattice position, there exists a bin. Once again, we 
label the bins by their serial indices k, k = 1,2,3, ... , M, 
placed from left to right. In the fe-th bin, there are a 
certain number of X particles, denoted by N Xfk and a 
certain number of Y particles, denoted by iV Yj fc. Both 
N Xy k and N Yi t are finite. The dynamics of the system is 
described by three basic processes: (i) Diffusion of the X 
particles in the fe-th bin to the (k— l)-th or the (k+ l)-th 
bin with a rate of diffusion 7. If an X particle in bin 1 
jumps towards the left, or an X particle in the M-th bin 
jumps to the right, then they are immediately replaced. 
(ii) Likewise, diffusion of the Y particles in the fc-th bin 
to the (fc — l)-th or the (k + l)-th bin with a rate of dif- 
fusion 7. If an X particle in bin 1 jumps towards the left, 
or an X particle in the M-th bin jumps to the right, then 
they are immediately replaced. (Hi) Reaction to produce 
an extra X particle having annihilated a Y particle (X + 
Y — > 2X), with a rate 1/N. To study the phenomenon 
of front propagation for this model, the initial configura- 
tion of the system is taken as iV Xj k = N[l — 0(fc — fco)] 
and JV Y , k = NQ(k — fc ) (a step function in the density 
of the X particles). Because of the reaction process (in), 
the number of X particles in any bin keeps increasing, 
until the supply of Y particles in that bin runs out. As a 
result, the region that is full of X particles slowly invades 
the region that is full of Y particles, and this constitutes 
a propagating front. 

The corres pon ding equation of the front that is analo- 
gous to Eq. (2.5), is slightly more complicated, and it is 
given by 



We now briefly turn our attention to the following 
reaction-diffusion process X+Y — > 2X on a lattice: at 



J 



£<JVx,k(t)> = 7 [(N x , k + 1 (t)) + (JV x , fc _i(t)) -2(N x . k (t)) 
at 



+ - [{N k (t)N x , k (t)) - «*(*)>] , (6-1) 



where N k (t) is the tot al n umber of particles in the fc-th 
bin at time t. In Eq. ( |0| ), if we replace (N k (t) N Xtk (t)) 
by (NM) (N Xtk (t)) = N(N Xik (t)), then one retrieves 
Eq. (2.5). In this section, therefore, our purpose is to in- 
vestigate if the correlation term (N k (t) N Xik (t)) has any 
bearing on the corrections of the asymptotic front speed 
over its corresponding value obtained from the model an- 
alyzed so far. 

Front propagation in this model has been studied nu- 
merically by Kessler and co-authors Go[. Our interest 
in this model is motivated by the following observation: 
in terms of the average number of X particles in a bin, 
an appropriate reaction rate yields an equation, which is 
similar to E q. (2.5). However, the linear growth term 
of Eq. (2.5) is replaced by a more complicated correla- 
tion function between the number of X and Y particles 
in the fc-th bin. Nevertheless, near the foremost bin of 
the X particles, the number of Y particles is so large that 
the fluctuations in their number remains small. Upon 



neglecting these fluctuations, the linear growth term for 
the X particles becomes the same as the one before, and 
one therefore expects to the speed correction to stay un- 
affected. Our purpose is to check this expectation nu- 
merically, by comparing data fo r th e front sp eed in this 
model with those given in Eq. (5.8) for Eq. (2.5). The 
algorithm that we use in our simulation is similar to that 
of Sec. V A. The value of M is taken to be 1500 and for 
the starting configuration of the system, we use ko = 50. 



The a symptotic front speed is calculated using Eqs. 
fvTi ">.7j. The measurement of the front speed starts at 
t = 200 and stops at t = 700. There are 985 Ai-values in 
this time interval. The corresponding v^(tj) vs. tj graph 
is shown below in Fig^ 1 1 . Using the same method of cal- 
culation as in Sec. VA, the asymptotic front speed for 
N = 10 and 7 = 1 comes out of the computer simulation 
as 



v N = 1.974 ±0.009. 



(6.2) 
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We notice that the error bars of Eqs. (5.8) and (p 



overlap with each other, and we conclude that the cor- 
relations between the total number of particles and the 
number of X particles in the bins [the (iVfc(t) N x , k(~t)) 
term in Eq. (|3.]|)] indeed do not affect the asymptotic 
front speed, as we had expected. 
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FIG. 11. Values of v N (tj) for 200 < tj < 700 and 
j = 1,2, ...,980. As expected, the fluctuations in iijv(tj) 
die out for large values of j. 



As noted before, Kessler et al. jl(J claimed that the 
prefactor of the speed correction was about a factor 2 
different from the value expected from the asymptotic 
formula of Brunet and Derrida. Whether this is 

due to the values of N not being large enough, or due 
to their particular way of implementing the stochastic 
simulations, is unclear to us. From our data, there is 
no reason to believe that for large N there is an essen- 
tial difference between the model with finite number of 
Y particles and the earlier model with an infinite supply 
of Y particles, otherwise the asymptotic formula would 
be incorrect as N — > oo. 



VII. CONCLUSION AND OUTLOOK 



In this paper we have identified a large number of ef- 
fects that play a role in the tip region of fluctuating 
fronts which in the mean-field limit reduce to pulled 
fronts. While a full theory from first principles, which 
yields explicit predictions for the front speed for finite N 
will be hard to come by, we believe that any such the- 
ory will incorporate most of the effects we have analyzed 
and studied with computer simulations. One important 
conclusion from our studies is that while the asymptotic 
large- N correction derived by Brunet and Derrida is uni- 
versal (in the sense of being independent of the details of 
the model) the corrections to this expression do depend 
on many details of the model. In most cases, deviations 



from the asymptotic results are significant for values of 
N that can realistically be studied. 

The message of this paper is as follows: the bulk of 
a fluctuating front can still be considered a uniformly 
translating one and one can properly define a co-moving 
co-ordinate system, in which the shape of the bulk re- 
mains unchanged; on the other hand, the position of 
the tip of the front in such a co-moving co-ordinate sys- 
tem fluctuates, and only on average the tip becomes sta- 
tionary in this co-moving co-ordinate system. From the 
mean- field limit of this fluctuating front, we know that 
the tip region is very important for its dynamics; as a 
result, the fluctuating tip plays a very significant role in 
deciding the asymptotic front speed, in which two very 
important aspects come to play a role - discrete nature 
of particles and discrete nature of the bin indices. In 
this paper, we have tried to formulate a theory to model 
this fluctuating tip. This theory is still a mean-field type 
theory. More specifically, at t = 0, the shape of the 
tip is always the same mean shape, and hence this the- 
ory is unable to capture the correlated diffusion events 
or the vacant foremost bin events (although we can es- 
timate the effect of the latter). Any alternative theory, 
that one might think at this point, must be able to take 
into account these fluctuation effects, which, as previ- 
ously explained, must be able to study multi-time corre- 
lation functions for correlated jumps at the tip. 

The prospect of such a theory however, looks grim at 
this point. Not only the problem becomes highly nonlin- 
ear, but also one must realize that the fluctuations in the 
number of X particles in the bins near the tip is of the 
same order as the number of X particles in them (~ 1), 
and there does not exist any small parameter that one 
can do perturbation theory with. 

Finally, we note here that we have confined our anal- 
ysis to cases in which the growth and hopping terms for 
few particles per site or bin are the same as those for a 
finite density of particles. In such cases the front speed 
converges to d* as N — ► oo. One should keep in mind, 
though, that there are also cases where with minor mod- 
ifications of the stochastic rules for few particles, one can 
arrive at a situation in which the speed does not converge 
to the pulled speed v* as N — ► oo, even though in the 
mean field limit one obtains a dynamical equation that 
admits pulled fronts. We will discuss this in more detail 
elsewhere p7j| . 

ACKNOWLEDGEMENT: D. P. wishes to acknowl- 
edge many important, motivating and stimulating discus- 
sions with Hans van Leeuwen, Ramses van Zon and Zohar 
Nussinov, the help from Ellak Somfai on setting up the 
computer simulations used to obtain results in this pa- 
per, and the financial support from "Fundamenteel On- 
derzoek der Materie" (FOM). 



23 



APPENDIX A: SUMMARY OF THE 
DERIVATION OF THE GENERALIZED 
VELOCITY CORRECTION FORMULA 

In this appendix, we der ive the generalized velocity 
correction formula, Eq. (3.4) and its interpretation. 

Without any loss of generality, we can express the front 
solution 0(0 for v N < v* by (Cf. §|]]|l||L7|]) 



0(0 = A sin [ Zi £ + [3 ] exp(- z r 



(Al) 



at the leading edge of the front, where z r — Re(z) and 
z.i = Im(z). The corresponding dispersion relation is then 
given by 



z r vn = 27 (cosh z r cos Zi — 1) + 1 
Zi vn = 27 sinh z r sin Zi . 



and 



(A2) 



The additive phase (3 in Eq. (3.3) can be scaled away by 
redefining A and the position of the origin from where 
£ is measured. We therefore drop f3 in this appendix. 
Since the scaled particle density has to be positive, i.e., 
0(0 > Oj the physical linear solution regime must be 
confined within the range where < z$ £ < tt. We now 
make a notational choice to denote the comoving coor- 
dinat e co rresponding to the node of the sine function in 
Eq. Q, by + 1, i.e., 



Zi (0 + 1) 



(A3) 



One should understand at this point that although Eq. 
(Al) suggests that there is a second node of 0(0; where 
the argument of the sin function becomes zero, such a 
node does not exists. Much before the argument of the sin 
function becomes zero, the nonlinear saturation term be- 
comes important and the solution (Al) for the linearized 
equation does not hold any longer. In this overly simpli- 
fied mean-field description, P^ys the role of the co- 
moving coordinate of the foremost bin. The mean field 
description of the front is then completed by claiming 
that 0(0 = for £ > + 1- Let, us also denote the 
density of the X particles in the "foremost bin" , which in 
this approximation is at 0, by a /N , to have 



A sin [ Zi ] exp(- z r 0) 



N 



(A4) 



Once the parameters A and a are known, Eqs. (A2 
A4) form a set of four equations for four unknowns, z r , 
Zi, an d vn, which we can then solve numerically for 
the asymptotic front speed vn- 

In order to put our results in a particular form that fa- 
cilitates comparison with th e earlier results in literature 
||, we analyze E qs (A2J-A4) for large N. First, with the 
help of the Eq. (A3), we reduce Eq. (A4) to 



A sin Zi exp(- z r 0) = 



(A5) 



Next, having introduced a new variable fi, such that 







In N 



ft. 



(A6) 



and using Eq. (A3), Eq. (A5) is further reduced to an 
implicit equation in /1: 



z r fl 



In 



.4 



In < sin 



In N + 1 + fi 



(A7) 



Since N is large, one can solve this implicit equation fi 
by means of a simple successive approximation proce- 
dure. At the lowest order, one can drop the /j, term in 
the d eno minator of the argument of the sine function in 
Eq. (A7) and obtain 



1 



In 4 

a 



In < sin ■ 



In N 



(A8) 



Finally, Zi can be obtained from Eqs. (A3), (A6) and 
©as 

(A9) 











7T Z r 


In N -) 


- Z r - 


M„ + 


{- In 








a 







By now, we have eliminated the unknown an d re- 
duced the problem to solving t hree unkn own s, z r , z^ and 
vn from three equations, Eqs. (A2) and (AS). From Eq. 
(|A9|), one can see that for large N, the approach of Zi to 
zero is extremely slow, going only as ln^ 1 N and also the 
fact that for the strict limi t of infinite N, Zi = 0, which 
reduces Eq. (A2) to Eq. (3.2), as it should. For large 
N, therefore, one expects that z r f=a zq and Zi ^ I, and 
one can expand vn around its value for z = zq . Upon ex- 
panding vn around v* , and using the solution of Zi from 
Eq. (|A9|) with z r replaced by Zq, we then find that the 



asymptotic speed vn is given by Eq. (3.4) 



APPENDIX B: DETERMINATION OF A 

Solving the full nonlinear difference-differential equa- 
tion, Eq. (2.6), is not an easy task. For a given set of 
values of 7 and N, there are essentially two methods to 
determine the value of A. The first one is to obtain the 
solution close to the saturation value 0=1 and thereafter 
iterate the solution until one reaches w 0. Close to the 
saturation value =1, one can define 77(0 = 1 — 0(0, 
which reduces Eq. (2J3) to an equation in 77(0, given by 



VN % = + + ^-!) - 2 ^)] 

-v(0 + v 2 (0- 



(Bl) 



For 77 values close to zero, the solution of Eq. (5/7) is 
given by the linearized equation 

- - TfoK + 1) + - 1) - 2r?(0] - r,(0 , (B2) 
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with the corresponding solution ry(£) = B \ exp[z;(£ — Co)]- 
Substitution of this solution in Eq. (B2) yields the dis- 
persion relation between vn and zi 



vnzi = 27 [ cosh zi — 1] — 1. 



One can then iterate this solution towards 



full solution of Eq. (Bl) can be written as 



B n exp[nz ( (£ - £0) 



(B3) 



0. The 



(B4) 



where the corresponding B n values are obtained from the 
recursion relation 



B,, 



B\B n -i + B 2 B n -2 + . . . + B n _iBi 
1 — nvNZi — 27 (cosh nzi — 1) 



(B5) 



As a starting point for constructing the solution near = 
1, one can choose arbitrary values of £0 an d B\, so long as 
the value of B\ is sufficiently small. At small values of B\ , 
any scaling of B\ amounts to a simple shift of the origin 
Co- Finally one can then match the solution, thus ob- 
tained, to the form 0(C) = Asin[zi(£ — £1)] exp[z r (£ — Ci)] 
near = and determine the value of A. 

The second method to obtain the numerical value of 
A is to assume a certain value of A close to <fi = with 
the functional form 0(C) = A sin[zi(£ — Ci)] ex P[ 2; r(C — Ci)] 
and then continue to iterate the corresponding solution 
in the direction of = 1 in a similar manner. This time, 
if the assumed value of A is correct, then close to = 1, 
one must recover the exponential behaviour of 0(C), as in 
Eq. (B4). However, we have found that the first method 
is stable under small changes in the starting value of 
Bi, while the second method is not stable under small 
changes in the assumed value of A. 

The first method should therefore be the natural 
choice, albeit from a practical point of view, one needs 
a very large number of B n values to extend the solu- 
tion of 0(C) all the way up to = 0. In practice, we 
have therefore used a "double shooting" method |p5| , in 
which the functions are calculated from both sides, and 
then matched somewhere in the middle. 

The matching of the values of the functions and their 
derivatives at 0o requires the values of zi, z r and Zi to 
be determined externally. For 7=1 and N = 10 4 , the 
values of z%, z r , Zi and A are numerically obtained from 
Eqs. |T~ 



(A.2) and (B3) as 



zi = 0.4187 T 0.0008, 
z r = 0.877 ± 0.002 , 
z t = 0.264 =p 0.007 and 
A = 0.961 ± 0.012. 



(B6) 



Of course, the numerical value of A depends on the ori- 
gin, where from C is measured for the form in Eq. (3.3). 
In Eq. (B6) above, the value of A is determined with 
/3 = 0, i.e., the solution of the linearized equation at the 



leading edge of the front is zero at £ — tt/ zi . We mention 
here that the uncertainty in these numbers is determined 
by the uncertainty in vjy, not by the inaccuracy of the 
numerical method. 
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The sequence of random numbers generated by drand48 
is random enough for our purposes so long as we use any 
large enough (~ 10 6 )initial seed. There is nothing spe- 
cial about choosing 123456, it is as arbitrarily chosen as 
any other initial seed ~ 10 6 . We have also run the sim- 
ulations for N = 10 4 and 7 = 1 with two other random 
seeds, namely s = 10 6 and s — 10 7 . For both of these two 



cases, the average speeds calculat ed by means of Eq. ( 5.7) 
fell within the error bars of Eq. (5.8). Thus, the chance 
of our results being affected by inherent correlations of 
the random number generator provided in the standard 
C library is ruled out. 
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